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ABSTRACT 

We present the results from the analysis of optical spectra of 31 Ha-selected 
regions in the extended UV (XUV) disks of M83 (NGC 5236) and NGC 4625 
recently discovered by GALEX. The spectra were obtained using IMACS at Las 
Campanas Observatory 6.5m Magellan I telescope and COSMIC at the Palomar 
200-inch telescope, respectively for M 83 and NGC 4625. The line ratios measured 
indicate nebular oxygen abundances (derived from the R23 parameter) of the 
order of Z q /5-Z q /1Cq For most emission- line regions analyzed the line fluxes 
and ratios measured are best reproduced by models of photoionization by single 
stars with masses in the range 20-40 M and oxygen abundances comparable 
to those derived from the R23 parameter. We find indications for a relatively 
high N/O abundance ratio in the XUV disk of M83. Although the metallicities 
derived imply that these are not the first stars formed in the XUV disks, such 
a level of enrichment could be reached in young spiral disks only 1 Gyr after 
these first stars would have formed. The amount of gas in the XUV disks allow 
maintaining the current level of star formation for at least a few Gyr. 

Subject headings: galaxies: abundances - galaxies: evolution - HII regions - 
techniques: spectroscopic - ultraviolet: galaxies 

1. Introduction 

Studying the outer edges of spiral galaxies provides a unique tool for understanding 
the formation and evolution of galactic disks. According to the inside-out scenario of disk 
formation, an increase in the gas-infall timescale with the galactocentric radius results in the 
delayed formation of the stars in the outer parts of the disk compared with its inner regions 
(Larson 1976, Matteucci & Francois 1989, Hou et al. 2000, Prantzos & Boissier 2000). This 
scenario, originally proposed to explain the color and metallicity gradients in the Milky Way, 
is also supported by recent N-body/SPH simulations of the evolution of galactic disks (e.g. 
Brook et al. 2006). 

Thus, star formation taking place now in some of the outermost regions of galaxies might 
probe physical conditions similar to those present during the formation of the first stars in 
the Universe. The low gas densities found in these regions also provide an excellent test 



1 We adopt 12+log(0/H) Q =8.69 (Allende Prieto, Lambert, & Asplund 2001) throughout the paper. 
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for the presence (or absence) of a threshold for star formation (Martin & Kennicutt 2001) 
and for studying the behavior of the star-formation law in low-density rarefied environments 
(Kennicutt 1989, Boissier et al. 2003, 2006, Elmegreen & Hunter 2006). 

Recent deep, wide-field observations of a sample of nearby spiral galaxies at UV wave- 
lengths carried out by the Galaxy Evolution Explorer (GALEX) satellite as part of its Nearby 
Galaxies Survey (Bianchi et al. 2003; see also Gil de Paz et al. 2006) have revealed the pres- 
ence of UV-bright complexes (XUV complexes hereafter) in the outermost parts of their 
disks. These galaxies, referred to as XUV-disk galaxies, host UV emission well beyond (~2- 
4x) their optical (D25) radii (Thilker et al. 2005a, 2007, in prep., Gil de Paz et al. 2005). To 
date, the two best-studied XUV disks are those of M83 (Thilker et al. 2005a; D=4.5Mpc) 
and NGC 4625 (Gil de Paz et al. 2005; D=9.5Mpc). These are two quite different exam- 
ples of XUV disks: M 83 has a very massive, large, high-surface-brightness optical disk with 
patchy XUV emission clearly disconnected from it; NGC 4625 is a low-luminosity system 
with ubiquitous XUV emission on top of an underlying low-surface-brightness optical disk 
(see Swaters & Balcells 2002). 

Any study of the properties of the XUV complexes discovered in M83 and NGC 4625 
must establish whether the UV emission associated with these extended disks is due to 
recent star formation or not. If that is the case, as previous studies have suggested, a 
detailed analysis of their properites should then provide fundamental clues to understand 
the possible mechanism(s) that led to formation of stars in these outermost regions of the 
disks. These can be internal mechanisms, such as density waves, or external influences, such 
as tidal interactions or satellite disruption. Finally, we should be able to establish whether 
(1) this is continously happening in the outer disks of these galaxies, (2) it is a transient but 
recurrent phenomenon that has taken place in these and other (perhaps all) spiral galaxies 
in the past, or, alternatively, (3) the XUV emission is a one-time phenomenon and these are 
the first generation of stars to form in the outer parts of these galaxies. 

In order to gain some deeper insight into the properties of these XUV complexes we 
have obtained deep optical spectroscopy of the outer disk of M 83 and NGC 4625 using the 
Magellan and Palomar 200-inch telescopes, respectively. In this work we study the physical 
conditions of the ionized gas (densities, metal abundances) in those XUV complexes showing 
Ha; emission. We also estimate the properties of the ionizing sources by comparing the line 
ratios measured with the predictions of photionization models; in particular, we determine if 
these line ratios are consistent with the line emission being due to photoionization by single 
stars as has been recently suggested (Gil de Paz et al. 2005). 

The spectroscopic observations of the XUV disks of M83 and NGC 4625 along with 
the data reduction methodology are described in Section [21 In that section we also briefly 
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describe complementary optical and GALEX UV imaging data. Section[3]presents the results 
from the analysis of the spectroscopic data. We discuss the nature of the XUV emission and 
its possible implications for the past and future evolution of these systems in Section |H The 
conclusions are summarized in Section [5J 

2. Observations and reduction 

2.1. GALEX UV imaging 

The regions analyzed here were originally identified in UV images taken by the GALEX 
satellite (Thilker et al. 2005a, Gil de Paz et al. 2005). The Galaxy Evolution Explorer 
(GALEX) satellite is a NASA's small explorer mission with a single instrument which consists 
of a Ritchey-Chretien telescope with a 50-cm aperture that allows simultaneous imaging in 
two UV bands/channels, FUV (A e fr=1516 A) and NUV (A e ff=2267 A), within a circular field 
of view of 1.2 degrees in diameter. See Martin et al. (2005) for a more detailed description 
of the mission and the GALEX instrument. 

The GALEX observations used in this paper were carried out on 2003 June 7 in the 
case of M83 and on 2004 April 5 for NGC 4625. The total exposure times (equal in both the 
FUV and NUV bands) were 1352 s and 1629 s for M83 and NGC 4625, respectively. These 
data were all reduced using the standard GALEX pipeline. 

2.2. Optical imaging 

Ground-based optical imaging has been carried out to obtain accurate positions for 
the regions responsible for the XUV emission. Accuracies of the order of a few tenths of 
arc second are needed for obtaining multi-object optical spectroscopy with relatively good 
spectral resolution for the 1-1.5 arcsec- wide slitlets used (see Section I2~3l) . 

On 2005 March 10-14, we obtained deep broad-band UR imaging data of the southern 
part of the XUV disk of M83 using the 2048x3150 pixels Site 3 Direct CCD at the Cassegrain 
focus of the Swope 1-m telescope at Las Campanas Observatory (Chile) (see Figure 1). The 
exposure time was 6 x 1800 seconds in U and 3x900 seconds in R. Deep ground-based optical 
imaging of NGC 4625 in B and R bands had been previously obtained at the Isaac Newton 
2.5-m telescope in La Palma (Spain) on 1995 May 28 and 1995 Dec 24 & 27 (see Figure 2). 
Images are available through the Isaac Newton Group (ING) data archive. See also Gil de 
Paz et al. (2005) and Swaters & Balcells (2002) for details on the optical imaging observations 
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of NGC 4625 (see also Table 1). 

Many of the galaxies with extended UV-bright disks found by GALEX so far show a 
relatively sharp cutoff in the azimuthally-averaged surface-brightness profile in Ha at the 
limit of the optical disk, while the UV light apparently extends smoothly beyond that radius 
(Meurer et al. 2004, Thilker et al. 2005a, Gil de Paz et al. 2005). It is therefore expected 
that a significant fraction of the regions identified in the GALEX images would not show 
any line emission and, consequently, could not be used for determining the ionized-gas metal 
abundances and dust reddening (from measurements of the Balmer-line decrement). With 
this in mind and in order to improve the efficiency of our spectroscopic observations, we also 
obtained deep narrow-band Ha imaging of both M83 and NGC 4625. M83 was observed 
on 2005 March 14 through a 70-A-wide narrow-band filter centered on 6570A using the 
2048x3150 pixels Site 3 Direct CCD attached to the Swope 1-m telescope in Las Campanas. 
The total exposure time in Ha was 10,200 seconds split in 9 exposures. On 2004 August 
20, we obtained an 800 second narrow-band image of NGC 4625 through a 20-A-wide filter 
centered on 6563A using COSMIC mounted at the prime focus of the Palomar Observatory 
5-m Hale telescope (Kells et al. 1998). 

Reduction of the images was carried out using For the Ha imaging data, the 

continuum was subtracted using a scaled i?-band image such that the fluxes of the field 
stars in both images were equal. The astrometry of the optical images was carried out 
using the IRAF tasks starfind, ccxymatch, ccmap, and ccsetwcs. We made use of 
the USNO-B1 catalog in the case of M83, and the USNO-A2 catalog for the images of 
NGC 4625. The r.m.s. achieved by our astrometric calibration was smaller than 0.3arcsec 
in all cases. The procedures followed for subtracting the continuum from the Ha images and 
for performing the astrometry of the optical images are described in detail in Gil de Paz, 
Madore, & Pevunova (2003). 

False-color RGB composites using GALEX FUV, broad-, and narrow-band ground- 
based optical images are shown in Figure 1 for M83 and in Figure 2 for NGC 4625. A 
summary of the UV and optical imaging observations is given in Table 1. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.3. Optical spectroscopy 

2.3.1. Sample selection 

Our spectroscopic sample is made up primarily of regions with detected emission in the 
continuum-subtracted narrow-band images. The vast majority of these regions correspond to 
stellar associations that also show bright UV emission. Additionally, we included [/-bright 
sources with no obvious Ha counterpart, mainly to take advantage of the large field-of- 
view offered by the IMACS spectroscopic observations of M83. In this paper we will focus 
exclusively on the analysis of the properties of the emission-line regions. 

In Table 2 we give the coordinates of the emission-line regions in the XUV disks of M 83 
and NGC 4625. These are all the regions for which our follow-up spectroscopic observations 
confirmed both their nature as emission-line sources (both Ha and H[3 being detected in 
emission) and having a redshift close to that of the corresponding parent galaxy. The co- 
ordinates given in this table are accurate to the level given by the r.m.s. of our astrometric 
calibration, i.e. better than 0.3 arcsec. Small systematic errors inherent to the USNO catalog 
could be also present. 

2.3.2. M83 observations 

Multi-object spectroscopic observations of the southern region of the XUV disk of M 83 
(see Figure 1) were obtained using the IMACS spectrograph at the 6.5-m Magellan-I (Baade) 
telescope in Las Campanas (Chile). The observations were carried out on 2005 April 18 using 
the f/4 long camera with the 300 1 mm -1 grating in its first order. This configuration provides 
a spatial scale of 0.11 arcsec pixel -1 , a reciprocal dispersion of 0.743A pixel -1 , and a spectral 
resolution of R~1200 for slitlets of 1 arcsec in width. IMACS uses a 8Kx8K-pixel mosaic 
camera. The objects observed were distributed over two different masks with total exposure 
times of 40 and 60 minutes. During the night the seeing ranged between 0.6 to 1.0 arcsec. 

The spectra obtained were reduced using the so-called COSMOS package, an IMACS 
data reduction pipeline developed by August Oemler. COSMOS consists of a set of tasks 
that allow correcting the science images for bias and flat-fielding, subtracting the sky, and 
performing a precise mapping of the CCD pixels onto the coordinate system of wavelength 
and slit position using the observing setup and mask generation files as input. 

In the case of the IMACS observations of M83 a total of 37 Ha-selected (34 in the 
XUV disk) and 32 [/-band-selected sources were included in the two masks observed. Out 
of the 37 Ha-selected sources 22 of them (19 out of 34 in the XUV disk) were detected 
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both in the Ha and H/3 lines in our spectroscopic data with redshifts close to that of M83. 
The positions and observed H/3 fluxes inside the slitlets for these 22 regions are given in 
Table 3. We also serendipitously discovered a background galaxy ([GMB2007]1) at z~0.3 
[RA(J2000)=13 h 37 m 24.53 s ; Dec(J2000)=-29°57'20.i"] whose emission in [OIII]A5007A falls 
in the window defined by the filter used in our narrow-band imaging observations and for 
which we also detected H/3 and [OIII]A4959A. 

2.3.3. NGC 4625 observations 

The spectroscopic observations of the XUV disk of NGC 4625 were carried out using 
COSMIC mounted at the prime focus of the Palomar Observatory 5-m (Hale) telescope. The 
spectra were taken on the night of 2005 March 16 using the 300 lmm' 1 grating and a mask 
with 1.5 arcsec- wide slitlets. COSMIC uses a single 2048 x 2048-pixels CCD detector. The 
total exposure time of the single mask observed was 40 minutes and the seeing during the 
observations was ~1. 2 arcsec. The spectra were reduced using standard IRAF tasks within 
the ccdred, onedspec, and twodspec packages. 

Out of the 17 Ha-selected targets included in the mask, 12 sources showed emission in at 
least Ha and H/3 at the redshift of NGC 4625 and were used in our analysis (see Table 3 for the 
list of positions and observed fluxes in H/3). Additionally, we identified a background galaxy 
QGMB2007]2) at ^-0.3 located at RA(J2000) = 12 h 41 m 41.89 s and Dec(J2000)=+41°16'10.9" . 

2.3.4- Flux calibration 

The spectroscopic observations of M 83 and NGC 4625 were both made under photomet- 
ric conditions. In the case of M83 we observed the spectroscopic standard star LTT 3218 
using a 2.5 arcsec- wide long-slit mask. We adopted the extinction curve for Cerro Tololo 
Inter-American Observatory. For the observations of NGC 4625 we obtained spectra of 
BD+33°2642 and HD 93521 using a 4arcsec-wide long-slit mask. The extinction curve of 
Kitt Peak National Observatory was adopted in this case. 

The flux-calibrated spectra obtained for the 22 (12) Ha-selected regions in M 83 (NGC 4625) 
used for this work are shown in Figure 3 (Figure 4). Since the observation of the spectropho- 
tometric standard stars was performed near the center of the field we compared the resulting 
flux-calibrated sky spectra obtained through each of the slitlets in order to ensure the validity 
of this calibration for all sources in the mask. 

For M83 the absolute calibration shows a dispersion of ~0.3mag (Figure 5a). However, 
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if we apply systematic offsets to the individual spectra (which do not affect the line ratios 
measured) the dispersion between these spectra due to changes in the wavelength dependence 
of the flux calibration is reduced to ~0.15mae; between 4000-6000 A and ~0.2 mag below 
4000 A and above 6000 A (Figure 5b). These numbers represent the average 1-a calibration 
errors for our spectra. 

In the case of NGC 4625 both the absolute and relative differences between all the sky 
spectra extracted are significanly smaller than for M83 (see Figures 5c and 5d). Once small 
global offsets are applied the r.m.s. between the spectra becomes ~0.1mag for A<4300 A or 
A>6600A or even smaller for 4300 A<A<6600 A. This is probably due to the much reduced 
field of view of this instrument and smaller format of the detector. 



3. Results 

3.1. Line fluxes and ratios 

The extinction-corrected H/5 line fluxes measured through our slitlets are given in Ta- 
ble 2. These spectroscopic fluxes all below the value expected for a region ionized by a single 
03 star (L Ha ^10 38 ergs" 1 ; Sternberg et al. 2003) at the distance of our targets. Similar 
conclusions were reached by Gil de Paz et al. (2005) through the analysis of narrow-band 
imaging data on NGC 4625. Note that although many of these regions are very compact, 
aperture effects consequence of the limited size of our slitlets (1-1.5 arcsec in width) might 
be significant. Therefore, the spectroscopic line fluxes given in Table 2 should be considered 
as lower limits to the total observed flux of the corresponding HII region. 

In Table 3 we present the line ratios measured for each of the 19 (12) confirmed line- 
emitting regions in the XUV disk of M83 (NGC 4625), plus three additional regions in its 
optical disk (NGC 5236: XUV 20, 21, and 22). Both fluxes and line ratios have been corrected 
for foreground Galactic and internal extinction using the Ha/H/9 Balmer-line decrement. 
We have adopted an intrinsic value for Ha/H/9 of 2.86 for case B recombination (Osterbrock 
1989) and the parametrization of the Galactic extinction law given by Cardelli, Clayton, & 
Mathis (1989). 

Figure 6 shows the standard diagnostic diagrams proposed by Baldwin, Phillips, & 
Terlevich (1981) and Veilleux & Osterbrock (1987). The line ratios measured in the XUV 
disk of M 83 (panels 6a & 6b) are similar to those found in local star- forming galaxies (SFG), 
represented here by the Ha-selected UCM-Survey galaxies of Zamorano et al. (1994, 1996). 
We find that regions in the inner (optical) disk of M83 (r<5.5'; filled circles) have lower 
excitations than those in the inner part of the XUV disk (5.5<r<10'; filled triangles) and 



- 9- 



even lower than in the outermost XUV disk (r>10'; filled stars). 

On other hand, the line-emitting regions in the XUV disk of NGC 4625 show relatively 
high excitation conditions. We also find a larger dispersion in the [OIII]A5007A/H/3 ratio for 
a given value of [NII]A6584A/Ha that in the case of M83. Gil de Paz et al. (2005) showed 
that many of the regions identified in the XUV disk of NGC 4625 have Ha luminosities 
consistent with photoionization due to single massive stars. Interestingly, simultaneous low 
[OIII]A5007A/H/5 and [NII]A6584A/H a ratios are predicted by some of the single-star pho- 
toionization models of Stasinska & Schaerer (1997). These models do not cover the entire 
parameter space of the properties of the ionizing star (luminosity, metallicity, etc.) or the 
surrounding HII region (density, geometry) nor were intended for reproducing the line ra- 
tios of single-star HII regions but their ionization structure. Thus, a detailed modelling of 
all the observed line ratios for each individual region covering a wider range of parameters 
than those explored by Stasinska & Schaerer (1997) is carried out in Section 13.31 (see also 
Figure 6). 

Panels 6b and 6d show the predictions of the photoionization models for evolving star- 
bursts developed by Stasinska & Leitherer (1996). In Figure 6b the models represent the 
variation in the line ratios with age between 1 to lOMyr for different metallicites (Z Q , Z /4, 
Z /1O, Zq/40). For starburst models with metallicities below solar the sequence in excitation 
conditions seen in these figures is mainly due to a change in age, from young high-excitation 
starbursts (whose effective temperature is dominated by very massive stars) at the top-left 
part of the diagram to more evolved low-excitation starbursts that can be found at the 
bottom-right. 

The fact that the loci of the emission-line XUV complexes in these diagrams coincide 
with those of local star-forming galaxies indicates that the line emission detected is produced 
by photoionization by young massive stars and not by planetary nebulae, which typically 
show lower [OII]AA3727AA/H/3 and somewhat higher [OIII]A5007A/H/5 line ratios than those 
found here (see e.g. Baldwin et al. 1981). The compact nature of these sources (see Figures 1 
& 2) also excludes the possibility of this line emission being associated with Diffuse Ionized 
Gas (DIG). 

3.2. Metal abundances based on strong lines 

Determining the metallicity of the ionized gas in star-forming galaxies and HII re- 
gions generally means measuring their oxygen abundances. Ideally this requires deriving 
the electron temperature from the measurement of a temperature-sensitive line ratio such 
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as [OIII]A4363A/[OIII]A5007A. Unfortunately, the [OIII]A4363A line is very weak and can 
only be detected in bright, actively star- forming low-metallicity systems. 

Since the [OIII]A4363A line in not detected in our spectra, in order to constrain the 
metallicities of the emission-line regions in the XUV disks of M83 and NGC 4625 we rely 
on spectral diagnostics based on strong lines such as those described in Kewley & Dopita 
(2002). Among these, perhaps the best estimator of the nebular oxygen abundance is the 
R23 ratio. Defined as log{([OII]AA3727AA+[OIIi]A5007A+[OIII]A4959A)/H / 3}, the R23 
parameter was first introduced by Pagel et al. (1979). Although very sensitive to metal 
abundance, this parameter is also sensitive to the ionization parameter, especially at low 
metallicities (Z<Z ). This degeneracy between metallicity and ionization parameter can 
be partially solved by taking into account the value of an indicator mostly sensitive to the 
ionization parameter, such as the [OII]AA3726,3729AA/[OIII]A5007A ([OII]/[OIII] hereafter) 
ratio (e.g. McGaugh 1991). 

Another critical issue regarding the use of the R23 ratio is the fact that is double- valued 
for most R23 flux ratios. It shows a maximum that, depending on the ionization parameter, 
is located anywhere between half solar to solar metallicity and below that maximum the value 
of the R23 can be reproduced by either a low-metallicity or a high-metallicity photoionization 
model. In order to solve this problem an initial guess on the oxygen abundance based on some 
other diagnostics or secondary indicators, such as the [NII]A6584A/[OII]AA3726,3729AA 
([Nil] /[Oil]), [NII]A6584A/[SII]AA6717,673lAA ([NII]/[SII]), or [NII]A6584A/H« ([NII]/Ha) 
ratios, is desirable. 

We have used two different recipes, those published by McGaugh (1991; M91 hereafter) 
and Pilyugin & Thuan (2005; PT05 hereafter), which allow for a determination of the oxygen 
abundance as a simple function of the R23 ratio and the [OII]/[OIII] line ratio, once either 
the low- or high-metallicity branch is chosen. Note that for high metallicities only the method 
of PT05 includes a dependence on the [OII]/[OIII] line ratio through the use of the excitation 
parameter, P{=([OIII]A5007A+[OIII]A4959A)/([OII]AA3727AA+[OIII]A5007A+[OIII]A4959A)}. 

3.2.1. M83 

The large dispersion in the line ratios measured across the disk of the M 83 indicates a 
wide range of ionization conditions and/or metal abundances. Therefore, we should be very 
cautious in choosing what branch of the parametrization of the R23 ratio to use. 

In the case of the three innermost regions analyzed in this paper (NGC 5236: XUV 20, 
21, and 22), which are still considered part of the optical disk of the galaxy, the [Nil] /[Oil] 
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ratios measured are high enough {log([NII]/[OII])>0 in all three cases} that we can confi- 
dently assume the high-metallicity-branch for the calibration of the R23 parameter. For the 
regions in the XUV disk of M83 the [Nil] /[Oil] ratios measured (between —1.1 and —0.3 in 
log scale) indicate metal abundances around solar. The [NII]/[SII] ratios measured, ranging 
between —0.2 and 0.4 in logarithmic scale but with most of the regions showing positive 
values, favors again metal abundances close to solar. In the case of the [Nil] /Ha line ratios, 
whose values range approximately between —1 and —0.4, we can only say that metallicities 
are probably above Z /2. We should point out that all [Nil] /[Oil], [NII]/[SII], and [Nil] /Ra 
are good diagnostics for the oxygen abundance as long as the N/O abundance ratio adopted 
by Kewley & Dopita (2002) is appropriate for these regions (see below for further discussion). 
These numbers do not enable us to determine which of the two branches of the calibration 
of the R23 ratio with metallicity should be used in the case of the emission-line regions of 
the XUV disk of M83. Accordingly, in Table 3 we give the values derived assuming both 
the low- and high-metallicity cases for both the M91 and PT05 calibrations. The fact that 
these secondary indicators suggest metal abundances close to the limit between the low- and 
high-metallicity branches of the R23 calibration results in very similar oxygen abundance 
estimates under the two assumptions for regions NGC 5236: XUV 07 and NGC 5236: XUV 
17. 

The R23 calibration of PT05 yields systematically lower oxygen abundances than those 
derived using the M91 method (similar systematic differences have been found by van Zee & 
Haynes 2006 when comparing the calibrations of M91 and Pilyugin 2000). This is the case for 
both the lower and upper branch of the calibration. These differences range between <0.1 dex 
in regions NGC 5236: XUV 01 and 09, NGC 4625: XUV 01 and 07 to >0.4dex in regions 
NGC 5236: XUV 07 and 17, and NGC 4625: XUV 06. The largest differences between the 
predictions of these two methods are found in regions located in the so-called transition zone 
(PT05). It is worth mentioning that the lower [upper] branch of the empirical calibration of 
PT05 was limited to oxygen abundances 12+log(O/H)<8.00 [12+log(0/H)>8.25] and values 
of the excitation parameter, 0.55<P<1 [P>0.1]. Thus, oxygen abundances in Table 3 not 
fulfilling these criteria are extrapolated values and should be considered with caution. 

The photoionization models of Kewley & Dopita (2002) assume a N/O abundance ratio 
that is approximately constant {log(N/0)~— 1.5; i.e. nitrogen is assumed of primary origin} 
at low oxygen abundances {12+log(0/H)<8.5} and rises linearly with the oxygen abundance 
above this value (i.e. secondary nitrogen). They used an empirical fit to the abundance 
measurements of van Zee, Haynes & Salzer (1998a). However, the N/O ratio is highly 
dependent on the actual stellar yields, star formation and gas-infall history of individual 
galaxies, and on the Initial Mass Function (IMF) (Prochaska et al. 2002; Pilyugin, Thuan, 
& Vilchez 2003; Koppen & Hensler 2005; Molla et al. 2006). In the case of the XUV disks 
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these properties may be different from those found in the HII regions of actively star-forming 
dwarfs and optical spiral disks, places where metal abundances are traditionally measured. 
For example, according to some of the scenarios proposed to explain the lack of line emission 
in XUV disks (Meurer et al. 2004), the IMF in these regions might have a lack of very massive 
stars. Such a top-light IMF would result in both the oxygen and nitrogen to be released by 
intermediate-mass stars which would led to a high N/O abundance ratio, higher than that 
adopted by Kewley & Dopita (2002) at low oxygen abundances. Should the actual N/O 
abundance ratio in XUV disks behave differently from what these authors assume some of 
the conclusions drawn above from the analysis of secondary indictors could change. A more 
extensive discussion on the possible N/O abundance ratios in these XUV disks is given in 
Section 13.3.21 

Moreover, some of the conclusions on the analysis of the secondary oxygen-abundance 
indicators are valid only if the ionization parameter in these regions is in the range considered 
by standard photoionization models such as those of Kewley & Dopita (2002). Here we should 
remind the reader that it has not been proved yet that photoionization models of typical 
HII regions adequately reproduce the physical conditions present in these, newly-discovered, 
extended UV disks. 

This paper is a first attempt to study the physical properties of the ionized gas in these 
regions, so we face many limitations that are inherent to previous empirical and theoretical 
studies mainly focused on the analysis of HII regions in the main body of spiral and irregular 
galaxies. In particular, the ionization parameter in the HII regions of the XUV disks could 
be either significantly higher or lower than in Galactic HII regions (due, for example, to a 
hotter/cooler mean stellar population). They could also show a different density structure 
than that found in typical HII regions and perhaps similar to that of the disks of Low 
Surface Brightness (LSB) galaxies, where the lower ISM pressure is thought to result in a 
more simple density structure (Mihos, Spaans, & McGaugh 1999). Moreover, single-star 
photoionization models might be more appropriate when trying to reproduce the conditions 
and output spectra of these regions (see Section |3~3]) . And finally, the N/O abundance ratio 
in XUV disks may be different from that measured elsewhere (see Section f3.3.2h . 

3.2.2. NGC 4625 

For NGC 4625 the [Nil] /[Oil] ratio in the emission-line regions of its XUV disk are 
all below or very close to —1.0 (in logarithmic scale), with the highest value being —0.94. 
Therefore, according to the photoionization models of Kewley & Dopita (2002), the metal- 
licities are necessarily below solar (for any ionization parameter between 5xl0 6 -or lower- 
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and 3x 10 s cms -1 ). We remind the reader that we are using 12+log(0/H)©=8.69 (Allende 
Prieto et al. 2001). Moreover, the 5-band luminosity-metallicity relationship recently de- 
rived by Salzer et al. (2005) indicates that for the luminosity of NGC 4625 (Mb—— 17. 4 mag) 
the vast majority of the local star-forming galaxies show oxygen abundances clearly below 
12+log(0/H)<8.5. 

Therefore, it is reasonable to conclude that the low-metallicity branch should be adopted 
for the XUV disk of NGC 4625, which is bluer and probably less chemically evolved than the 
overall galaxy stellar population (Gil de Paz et al. 2005). It is worth emphasizing here that 
the [Nil] /[Oil] ratio depends only weakly on the ionization parameter. This conclusion would 
still be valid even if the N/O abundance ratio in these regions is higher than log(N/0)=— 1.5, 
the constant value adopted by Kewley & Dopita (2002) for primary-nitrogen production. 

In one of the regions of the XUV disk of NGC 4625 (NGC 4625: XUV 09), the R23 
value derived is somewhat larger than the maximum value of this parameter for the cor- 
responding [OII]/[OIII] ratio (or excitation parameter, P), which results in an abnormally 
high value for the oxygen abundance, in spite its low [Nil] /[Oil] ratio (—1.2 in the log). 
This suggests that in this region the metal abundance is probably close to the value for 
which the R23 ratio is maximuum according to both photoionization models and empirical 
recipes {12+log(0/H)~8.4 for log([OII]/[OIII])~0; see e.g. Kobulnicky et al. 1999} and that 
observational errors might have moved its value above this maximuum. Therefore, in the 
particular case of NGC 4625: XUV 09 we adopt 12+log(0/H)=8.4 as our best metallicity 
estimate. The PT05 calibration is of limited use in the case of NGC 4625, since all regions 
analyzed in its XUV disk (except for NGC 4625: XUV 01) show excitation parameters, P, 
well below 0.55, the minimum value given by PT05 as giving reliable results. Its PT05-based 
oxygen abundances are, consequently, extrapolated values. 



3.3. Photoionization modeling 

3.3.1. Physical parameters of the models 

In order to study the physical conditions in the XUV disks we have compared the 
observed line ratios with the predictions of a complete set of ~13,000 photoionization models 
covering a wide range of physical properties both for the gas and for the source of ionizing 
radiation. Calculations were performed with version 06.02 of CLOUDY, last described by 
Ferland et al. (1998). 

First, we ran two sets of 3,262 simulations each using single massive stars from the 
CoStar library of Schaerer & de Koter (1997) with ZAMS masses 20, 25, 40, 60, 85, and 
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120 M , ages between 1 and 7Myr (or the lifetime of the star, whichever happens first) 
with a step of 1 Myr for stars >25 M and of 2 Myr below that mass, and two different 
stellar metallicities (Z and Z /5). The luminosity of the star was determined using the 
relations given by Tout et al. (1996). For the gas we adopted a spherical distribution with 
constant density 1, 10, 10 2 , and 10 3 cm -3 and a inner radius between 10 16 cm to 10 17 5 cm (in 
steps of 0.5 dex). The outer radius was set to 10 21 cm (~325pc) in order to ensure that all 
modeled HII regions were ionization bounded. The metal abundance of both gas and dust 
grains was allowed to simultaneously vary between Z /4O and 2xZ , with values (1/40, 
1/20, 1/10, 1/5, 3/10, 2/5, 1/2, 3/5, 4/5, 1, 1.5, and 2xZ ). The first set of realizations 
was run adopting a constant solar N/O abundance ratio {log(N/O) =— 0.76} while for the 
second one we adopted a constant value of log(N/0)=— 1.5 for oxygen abundances below 
12+log(0/H)=8.45 and a N/O rising linearly with O/H above that value, i.e. a behavior 
similar to that adopted in the photoionization models developed by Kewley & Dopita (2002). 
We stopped the calculation when the gas temperature reached 3,000 K, value at which the 
contribution of the gas to the optical line emission is negligible. 

As we commented in Section the weak line fluxes of these HII regions indicate that they 
are most likely powered by single, at a maximum, a few massive stars and not by a massive 
stellar cluster. In order to further verify this result we will also ran two sets of models (with 
two different N/O abundance ratios) of 3,456 CLOUDY realizations each where the gas was 
ionized by a stellar cluster with a mass of 10 3 M (see Gil de Paz et al. 2005). For the stellar 
continuum we used the predictions of the Starburst99 population systhesis models (Leitherer 
et al. 1999) for metallicities between Z /2O and 2xZ (with values 1/20, 1/5, 2/5, 1, 2xZ ; 
i.e. those of the Geneva tracks used in these realizations) and ages in the range 1-9 Myr (with 
steps of 0.1 Myr), assuming the IMF of Kroupa, Tout, & Gilmore (1993). We made use of 
the Starburst99 default (recommended) stellar atmospheres of Pauldrach/Hillier (Pauldrach 
2005, Hillier 2003, Lejeune, Cuisinier, & Buser 1997). The spatial distribution and chemical 
composition of the gas and dust grains were those used for the single-star photoionization 
realizations described above. We used a somewhat finer grid of gas densities (0.5 dex steps 
instead of 1 dex) in order to better explore the ionization parameter. Thanks to the homology 
relations of photoionization models with the ionization parameter the use of a fixed stellar 
mass for the cluster should not strongly limit the range of possible line ratios derived as long 
as the density of the gas is adequately sampled. 

For ~10% of the starburst models the outer radius (10 21 cm, i.e. 325 pc) was reached be- 
fore the temperature was below 3,000 K. This happens only for very low densities (io ~ - 5 cm~ 3 ) 
and only for low gas metallicities or for high metallicities when the ionizing starburst is young 
(<3Myr). None of the best-fitting solutions derived correspond to one of these "density- 
bounded" models. 
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The best-fitting set of models in each case was determined by ^-minimization of the 
extinction-free line intensities normalized to H/? in logarithmic scale. We adopted an average 
error of 0.10 dex in all line ratios, except for [OII]AA3726,3729AA/H/3 and [SII]AA6717,673lAA/FT3 
where errors of 0.15 and 0.20 dex were assumed, respectively (see Section [2.3.41) . We checked 
that the intrinsic Ha/H/? ratio predicted by these realizations was similar (within ~5%) to the 
case-B recombination coefficient assumed as part of our procedure for correcting the observed 
line ratios for extinction. Only the properties of regions for which the [OII]AA3726,3729AA/H/3, 
[OIII]A5007A/H/?, and [NII]A6584A/H/3 line ratios are available were determined. When the 
[SII]AA6717,673lAA/H^ line ratio was available, we included it also in the minimization. We 
note that the S/O abundance ratio is approximately constant with the oxygen abundance 
and similar to the solar value (see Kehrig et al. 2006 and references therein). 

3.3.2. N/O abundance ratio 

The N/O abundance ratio is a function of the stellar yields, the star formation and 
gas-infall histories of galaxies, and of the IMF. Most galaxies follow a relatively well defined 
sequence in the N/O-O/H diagram with a constant value of the N/O ratio for oxygen abun- 
dances below 12+log(0/H)=8.5, which suggests a primary origen for the nitrogen, and a 
N/O increasing linearly with O/H above that value, which is explained in terms of a sec- 
ondary production of nitrogen (van Zee et al. 1998a; Liang et al. 2006; see also Molla et al. 
2006 for a theoretical review on the subject). However, as we commented in Section [372TTJ 
both the history and IMF of the outermost regions of spiral disks (such as those analyzed 
here) may be quite different from those of the galaxies where nitrogen and oxygen abun- 
dances have been traditionally measured. Due to this limitation we decided to consider two 
possible N/O abundance ratios in our photoionization models: one similar to that used by 
Kewley & Dopita (2002) and that is based on the results of van Zee et al. (1998a) and other 
where N/O is constant throughout the entire range in oxygen abundance and equal to the 
solar value. The main difference between the predictions of these two sets of models is the 
intensity of the [NII]A6584A line. This is due to the fact that this line or other nitrogen lines 
have only a small contribution to the cooling of the gas, so a change in the N/O abundance 
ratio does not alter significantly the temperature structure of the nebula. 

Tables 4 through 7 show that in the case of M 83 the models where the N/O is assumed 
to be constant and equal to the solar value provide a better agreement to the data than 
those obtained using a O/H-scaled N/O abundance ratio. The observed increase in the 
best-fitting reduced x 2 f° r M83 is consequence of the [NII]A6584A line in the models being 
too weak for any set of physical parameters of the nebula compatible with the observed 
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oxygen line ratios. This can be clearly seen in Table 8 where we compare the observed line 
ratios (in logarithmic scale) with those predicted by the best-fitting photoionization model 
obtained for each combination of ionization source and N/O ratio. The models with O/H- 
scaled N/O abundance ratio systematically underpredict the value of the extinction-corrected 
[NII]A6584A/H/? line ratio. 

This result suggests that the actual N/O abundance ratio in the XUV disk of M 83 might 
be higher than that obtained in actively star-forming dwarf galaxies and in the optical disks 
of spiral galaxies at low oxygen abundances {log(N/0)=— 1.5}. There are several possible 
explanations for such a high N/O, e.g. the presence of a dominant evolved stellar population 
(Pilyugin et al. 2003; Molla et al. 2006) or a top-light IMF. The former scenario seems 
unlikely considering the relatively blue colors measured in the outer parts of these and other 
spiral disks (Munoz-Mateos et al. 2006; see also Section fO|) . Alternatively, an episodic star 
formation history with the last episode of star formation taking place several Gyr ago (i.e. 
time long enough for the release of nitrogen by intermediate-mass stars) could also result in 
relatively high N/O ratios. 

Another possible scenario involves the infall of pristine gas in regions with moderate 
(close to solar) oxygen abundance where the nitrogen production is already secondary. This 
gas infall would drive the system into a region of lower oxygen abundance without changing 
the N/O ratio (Koppen & Hensler 2005). This scenario has been also proposed to explain 
the large dispersion in the N/O ratio at oxygen abundances 12+log(0/H)>7.9 (including 
some objects with high values of N/O) recently found in a sample of dwarf irregular galaxies 
with low ionization parameters (van Zee & Haynes 2006). 

In the case of NGC 4625 the results obtained for the two different N/ O ratios considered 
are comparable except for region NGC 4625: XUV 05 where the solar N/O yields a better 
fit and for regions NGC 4625: XUV 09 and 12 where a O/H-scaled N/O ratio is favored. 
However, for most regions (all except NGC 4625: XUV 01) the line fluxes predicted by the 
best-fitting single-star models with a O/H-scaled N/O ratio are systematically lower than 
the lower limits measured by up to a factor of 40 (see Tables 2 & 5). 

For the remaining of the paper we will focus on the results obtained for the set of models 
with solar N/O abundance ratio (Tables 4 & 6) since they provide a better fit to the data 
overall. 
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3.3.3. Best-fitting single-star models 

In Table 4 we show the parameters of the best-fitting single-star photoionization model 
with solar N/O ratio. We find that the lowest reduced x 2 achieved are not very different 
from 1 for most of the regions, which indicates that these single-star photoionization models 
provide a good fit to the line ratios measured. The fact that they are not much smaller 
than 1 also indicates that we are not over-fitting our data. The high reduced \ 2 found in 
some of the regions might be due to these regions being powered by a few stars of different 
effective temperature or by a larger composite stellar population such a stellar cluster (see 
below). We have estimated the uncertainties in the properties derived by determining the 
minimum (subscripts in Tables 4 & 5) and maximum (superscripts in Tables 4 & 5) values 
of the each parameter that result in values of the absolute (i.e. not reduced) x 2 that difer 
from the minimum \ 2 by less than A% 2 =1 (see Avni 1976). This yields the 68.3% confidence 
intervals in each parameter when they are considered separately. Table 4 shows that while 
the ages and masses of the ionizing stars and the metal abundance of the gas are relatively 
well constrained, the inner radius of the gas distribution is very poorly constrained by our 
measurements. 

In roughly half of the regions (NGC 5236: XUV 01, 07, 10, 11, 17, and NGC 4625: 
XUV 01, 05, 09) the H/3 fluxes predicted by the best-fitting model are comparable or a 
factor of a few above (explicable by slit-aperture effects) the fluxes measured through the 
slits (Table 2). There are some cases (NGC 5236: XUV 09, 12, 14, and NGC 4625: XUV 02, 
06, 07) where the measurements are above the models but only by a factor ~2-5x. These 
HII regions could well be powered by a very few massive stars. Finally, in a few regions 
(NGC 5236: XUV 13, 20, 21, 22, and NGC 4625: XUV 06, 12) the slit R(3 fluxes are more 
than a factor of 10 above those predicted by the models suggesting that the source of the 
ionizing radiation is a relatively massive stellar cluster. This was expected in the case of 
regions NGC 5236: XUV 20, 21, and 22 since they are very luminous regions in the optical 
disk of M83. 

In the case of M 83 the ZAMS masses derived for the ionizing stars of the XUV regions 
are in the range 20-40 M with a best-fitting age of 1 to 3 Myr. With respect to the properties 
of the gas, the densities and inner radii show a wide range of values between 1-10 3 cm -3 and 
10 16 -10 17 5 cm, respectively. The best-fitting single-star models for NGC 4625 yield lower gas 
densities and somewhat higher temperatures and masses for the ionizing stars than those of 
M83. 

The gas metal abundances nicely agree with those derived in Section 13.21 using the R23 
ratio. Indeed, while M83 shows a strong metallicity gradient with abundances between 
Zq/10 and Z /4 in the outer disk and close or even above solar in the optical disk, the 
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XUV-disk regions of NGC 4625 (except for NGC 4625: XUV 09, with x 2 >5) show best- 
fitting metallicities that are either Z /5 or Z /1O. 

In light of these results we can now confidently say that metal abundances as low 
12+ log (0/H)=7.86 (regions NGC 5236: XUV 01 and NGC 5236: XUV 10) and 12+ log (0/H)=7.94 
(NGC 4625: XUV 01) are found in the XUV disks of M83 and NGC 4625, respectively. 
These metallicities are comparable to the lowest oxygen abundances found by van Zee et 
al. (1998b; 12+ log (0/H)=7.92 in NGC 5457) and Ferguson, Gallagher, & Wyse (1998; 
12+ log (0/H)=7.95 in NGC 1058) in the outermost regions of a sample of spiral galaxies. 
Note that in the case of Ferguson et al. (1998) only galaxies in Ferguson (1997) that showed 
very extended HI and extreme outer-disk star formation were analyzed. 

3.3.4- Best-fitting starburst models 

Table 6 shows the properties and reduced x 2 °f the best-fitting starburst models with 
solar N/O ratio. The reduced \ 2 values are larger (i.e. worse than) or, at best, comparable 
to those obtained from the corresponding single-star photoionization model. This indicates 
that for low stellar masses such as those of the XUV complexes (Gil de Paz et al. 2005) 
the shape of the ionizing spectrum is different from that of a modeled composite stellar 
population, which are usually built for massive star-forming regions and starburst galaxies. 
In particular, even in those cases where the presence of more than one ionizing star is 
suspected, stochastic effects on the upper IMF might result in the ionizing radiation from 
these regions to be dominated by only two or three of these stars, which would lead to an 
ionizing spectrum different from that of a low-mass stellar cluster with a continuous IMF, 
such as those generated by the Starburst99 population synthesis models. Unfortunately, 
the constraints imposed by the few line ratios measured in this study are limited so the 
determination of the number and masses of the ionizing stars present in each of these regions 
is beyond our reach. High-resolution imaging observations using HST would partially solve 
this problem by resolving individual massive stars at the distance of the nearest XUV disks. 
As expected, younger best-fitting starburst models typically correspond to regions where the 
effective temperature of the best-fitting single-star model is higher. The 68.3% confidence 
intervals in metal abundance of some of the regions in NGC 4625 are noticeably wider that 
those found using single-star models. Again, we find that the inner radius of the dust 
distribution is very poorly constrained. 

We note that although for the low metal abundances of the regions analyzed here the 
differences between the CoStar library and Pauldrach/Hillier stellar atmospheres used in 
Starburst99 are not as dramatic as for Z>Z (Smith, Norris, & Crowther 2002) differ- 
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ences between the two models could also partially explain the different best-fitting \ 2 values 
achieved. Should that be the case the massive stars present in the emission-line regions in 
the XUV disk would likely have a harder ionizing spectrum (as in the CoStar library; see 
Smith et al. 2002) than the O-stars included in the Starburst99 models. 

However, we should emphasize here that the fact that single-star models better repro- 
duce the line ratios was suspected already based uniquely on the low emission-line fluxes 
measured in these regions (see Table 2). Unfortunately, as we commented in Section f3.3.1[ 
some of these spectroscopic fluxes might be significantly affected by aperture effects, so our 
results could not be based solely on that fact. 

3.4. Metallicity gradients 

The oxygen metal abundances of the HII regions in the disk of M 83 are found to decrease 
with the galactocentric distance (Figure 7a). If we assume that for regions NGC 5236: 
XUV 09 to NGC 5236: XUV 14 the low-metallicity branch of the R23 calibration is the 
most appropiate the best linear fit yields a metallicity gradient of d[0/H]=— 0.112 dexkpc -1 
(d[O/H]=-0.069dexkpc- 1 ) using the M91 (PT05) method, with 

12 + log (O/H) = -0.112 x r(kpc) + 9.36; a = 0.37dex (M91) (1) 
12 + log (O/H) = -0.069 x r(kpc) + 8.60; a = 0.34 dex (PT05) (2) 

Here we have combined our metallicity measurements for the XUV disk with those derived 
from the line ratios of HII regions in the optical disk (Webster & Smith 1983) using the 
same methodology as for the XUV-disk abundances (open dots in Figure 7a). If the high- 
metallicity branch of the R23 calibration is now assumed for regions NGC 5236: XUV 09 to 
NGC 5236: XUV 14 the metallicity gradient is much shallower, 

12 + log (O/H) = -0.051 x r(kpc) +9.26; a = 0.16 dex (M91) (3) 
12 + log (O/H) = -0.017 xr(kpc) + 8.51; a = 0.25 dex (PT05). (4) 

Finally, if the branch of the R23 calibration is assigned based on the metal abundance of the 
best-fitting CLOUDY single-star photoionization model with solar N/O, i.e. low-metallicity 
branch for all regions except for NGC 5236: XUV 12, 14, 20, 21, and 22, the result is 

12 + log (O/H) = -0.100 x r(kpc) +9.41; a = 0.30 dex (M91) (5) 
12 + log (O/H) = -0.058 xr(kpc) + 8.64; a = 0.30 dex (PT05). (6) 

These values are slightly stepper and shallower, respectively, than the metallicity gradient 
measured for the Milky Way (d[0/H] M w~~ 0.08 dexkpc -1 ; Boissier et al. 1999 and references 
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therein) and than those derived for a sample of spiral galaxies with oxygen abundances 
measurements at large radii (van Zee et al. 1998b, Ferguson et al. 1998). If we now express 
the metallicity gradient in units of the blue-band scale- length (r^) or D25 radius (r 2 s) we 
obtain —0.198 dex/r^ and —1.01 dex/r 2 5, respectively, in the case of the M91 calibration and 
—0.115 dex/Vd and — 0.59dex/r25 with the PT05 calibration. These values are similar to 
the "universal" gradient of — 0.2dex/r<2 (or — 0.8dex/r 2 5) proposed by Prantzos & Boissier 
(2000, see also Garnett et al. 1997) for non-barred disk galaxies and within the range of 
values found by van Zee et al. (1998) and Ferguson et al. (1998) in spirals with outer-disk 
star formation. 

In light of Figure 7a instead of a single linear fit, the radial dependence of metal abun- 
dance seems to suggest a smooth metallicity gradient below lOkpc and sharp decrease in 
metallicity from ~Z© to Z /5 (if the M91 calibration is used) at that radius. This behavior 
could be explained if there exists a cut-off radius beyond which the disk of M 83 is signif- 
icantly less efficient in forming stars and it is allowed to evolve mostly viscously (Clarke 
1989). Note, however, that this might also be the consequence of the large uncertainties 
expected for any R23-based metallicity estimate at abundances close to half the solar value. 
A larger number of metallicity measurements, specially in the range 10-15 kpc, and more 
precise estimates of the metallicity using the [OIII]A4363A temperature-sensitive line would 
be required to confirm (or rule out) this. 

In the case of NGC 4625, where the low-metallicity branch can be safely assumed, the 
best linear fits obtained are 



The metallicity gradient derived, A[0/H]/Ai?=— 0.037±0.032 dexkpc -1 (using the M91 cal- 
ibration), is very modest. Even after normalizing this value by the galaxy scale-length or 
D25 radius, this is one of the shallowest gradients measured in disk galaxies, —0.080 dex/r 2 5. 
Note, however, that NGC 4625 is significantly less luminous (M B =— 17.4 mag) than the spiral 
galaxies where metallicity gradients have been traditionally measured. 

The dust attenuation also seems to decrease with the galactocentric distance in the 
case of M83 although with a very large dispersion (Figure 7b). The fact that this scatter 
is present even when neighboring regions are considered (e.g. NGC 5236: XUV 07 and 
NGC 5236: XUV 08) suggests a patchy dust distribution. In the case of NGC 4625 (see 
Figure 7d) no obvious trend with the galactocentric distance is found. 

It is also remarkable that positive and relatively large values of C(H/3) are measured 
very far out into the disks, showing the presence of dust in these galaxies at distances larger 



12 + log(0/H) 
12 + log (O/H) = 



-0.037 x r(kpc) + 8.44; a = 0.20 dex (M91) 
0.017 x r(kpc) + 8.16; a = 0.39 dex (PT05). 



(7) 
(8) 
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than 18kpc in the case of M83 and more than lOkpc in NGC 4625 (see Popescu & Tuffs 
2003 for an example of dust emission associated with an extended HI disk). 

4. Discussion 

4.1. Nature of the XUV disks: XUV complexes 

From the analysis carried out in this paper and from previous studies on the extended- 
UV disk phenomenon (Thilker et al. 2005a, Gil de Paz et al. 2005) is now clear that the 
XUV emission is due to young stars associated with low-mass stellar associations located at 
large galactocentric distances. The line ratios measured in those XUV complexes showing 
line emission indicate the presence of UV ionizing radiation being emitted locally in these 
regions. It is therefore unlikely that the non-ionizing UV light responsible for the XUV 
emission seen by GALEX is due to scattering by dust, like it is the case of the polar regions 
of the starburst galaxy M 82 (Hoopes et al. 2005). Our measurements also exclude planetary 
nebulae (with very high temperatures and highly-excited gas; see Section 13.3.11) as being the 
agents responsible for the observed XUV emission. Finally, a significant contribution by blue- 
HB stars it is highly unlikely considering (1) the blue NUV— optical colors of the XUV disks, 
(1) the high equivalent widths measured in some XUV complexes and (2) the morphology 
(resembling spiral arms) of the XUV emission; all unusual circumstances in objects with UV 
emission dominated by evolved stars (globular clusters, elliptical galaxies) but common in 
galaxies forming massive stars in recent epoch. 

Our results also support the conclusions of Gil de Paz et al. (2005) that suggested that 
the majority of the admittedly scarce emission-line sources found in the XUV disks of these 
galaxies seem to be powered by single stars, with masses between 20-40 M Q in most cases. 
As a consequence of this, a stochastic treatment of the IMF (e.g. Cervino 1998) is necessary 
in order to properly analyze the properties of the stellar populations in these XUV disks and 
perhaps in the outer edges of all spiral galaxies. In particular, knowledge on the fraction of 
low-mass stellar complexes that, once stochastic effects are properly accounted for, would 
form at least one ionizing star is key to understanding the observed differences between the 
UV and Ha light profiles of XUV disks (see e.g. Meurer et al. 2004). 

In this sense, it has been recently argued (Boissier et al. 2006) that the lack of Ha 
emission in the azimuthally-averaged profiles of spiral galaxies beyond a given radius, which 
is commonly associated with the existence of a intrinsic star formation threshold (e.g. Martin 
& Kennicutt 2001), could instead be a consequence of the fact that at those low levels of 
star formation the number of ionizing stars expected to be found (at a given time) inside 
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one of these annuli of the radial profile is less than unity. Under these circumstances, the 
chances of catching one of these ionizing stars are very low, leading to an effective lack of 
HII regions in most cases, but not necessarily to a lack of time-integrated star formation. 

4.2. Nature of the XUV disks: Extended HI, cause or consequence? 

Another property of the XUV disks is that they are always found in galaxies with very 
extended HI distributions (see Tilanus & Allen 1993 in the case of M 83 and Bush & Wilcots 
2004 for NGC 4625). This is not surprising since gas is of course a necessary ingredient 
for the formation of the stars responsible of the UV emission and the neutral-to-molecular 
gas ratio is thought to increase with the galactocentric radius in field spiral galaxies (e.g. 
Kohta, Naomasa, & Kuno 2001). However, the fact that significant amounts of HI (compared 
with their stellar content) are still found at the positions of these UV complexes suggests 
that the efficiency in converting neutral gas into stars is relatively low and that, under the 
right circumstances, these XUV disks could maintain the observed level of star formation 
for a cosmologically significant period of time. The SFR per unit area of these XUV disks 
can be obtained from the azimuthally-averaged FUV-surface-brightness after a correction 
for extinction of A F uv=l ma g is adopted^. If we now use the gas- surface- density profile of 
Thilker et al. (2005a) the consumption timescale derived for M83 is >6 Gyr at all radii, only 
considering the HI. In the case of NGC 4625 the consumption timescale is ~3 Gyr in the 
innermost parts of the XUV disk and rises up to values of the order of a Hubble time in the 
outer regions (Gil de Paz et al. 2005). 

Alternatively, it could be argued that the presence of an extended HI disk is a con- 
sequence of the XUV emission instead of a necessary pre-existing condition (Allen 2002). 
Under this scenario, the extended HI emission would be the result of the dissociation of H2 
in photodissociation regions (PDRs) by the UV photons emitted by the same newly-formed 
stars responsible for the XUV emission. 

However, in most of these XUV disks the complexes responsible for the FUV emission 
cover only a fraction of the area of the HI disk. Therefore, in order for this scenario to be valid 
in the case of the XUV disks the time for the reformation of the H2 molecule on the surface 
of dust grains from atomic hydrogen, Tf orm , should be larger than, or at least comparable 
to, the average time a particular region in the disk is embedded in FUV radiation. This 



3 This corresponds to an average observed (FUV— NUV) color of 0.3 mag and the relationship between 
(FUV— NUV) color and Apuv given by Boissier et al. (2006). Note that UV magnitudes are in AB scale and 
optical magnitudes are referred to Vega throughout the text. 
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condition can be written as 

^ ^HI / n x 

Tform £ X T FUV (9) 

^iFUV 

where Am is the area of the XUV disk where HI emission is detected, Apuv is the area of 
the disk reached by the FUV photons (responsible for the H 2 photodissociation in PDRs), 
and tfuv is the timescale of the FUV emission, which is the order of a few hundred Myr for 
an instantaneous burst. 

The H_2 reformation rate can be obtained using the expression given by Hollenbach & 
McKee (1979), 

Rh 2 = = ^n u v U (T^r]n g (10) 

where tih is the number density of hydrogen atoms (in atomic form), v# is their themal 
speed, a is the cross section of the grain, £ is the sticking coefficient, i] is the probability of 
bond formation once two hydrogen atoms encounter each other and n g is the number density 
of dust grains in the ISM. Here we assume that the contribution of the H~ process to the 
formation of H 2 is negligible. From the expression of the reformation rate given above, the 
reformation time, defined as the time necessary for half of the total hydrogen atoms to be 
in molecular form, can be written as 

i form 7- v / 

For simplicity, we assume a gas temperature of 70 K (Hughes, Thompson, & Colvin 1971), 
which yields a thermal speed for the hydrogen atoms of 1.2xl0 5 cms -1 (Spitzer 1978). The 
cross section of the grain is computed as 7ra 2 , where a is the typical radius of the grains 
in the ISM (0.1 /im; Jones & Merrill 1976). The sticking coefficient is adopted to be £=0.3 
and Tj—1 (Vidali et al. 2005; Hollenbach & Salpeter 1971). The number density of grains is 
estimated to be 210~ 12 x-n,H after an average grain material density of 2 gem -3 (Draine & 
Lee 1984) and the Galactic ISM gas-to-dust mass ratio are adopted (i? 9 /d~100; Knapp & 
Kerr 1974). This yields a reformation time of 

10 3 

r /orm (Myr) « (12) 

n H 

where nn is expressed in cm -3 . Allen (2002) gives a similar expression although his refor- 
mation timescales are somewhat smaller than those given by Equation [12] This is in part 
due to the fact that he uses parameters typical of PDRs while in our case the reformation 
of the H 2 is expected to take place in a more rarefied ISM. 

For XUV disks where the FUV emission covers a high fraction of the disk (like in the 
case of the inner XUV disk of NGC 4625; Gil de Paz et al. 2005) the reformation time might 
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fulfill the condition imposed by Equation [9] as long as the gas is not very dense. On the 
other hand, in M83, where the FUV emission covers only ~10% of the extended HI disk, the 
gas density required for the hydrogen to stay in atomic form is very low (nn<l cm" 3 ). We 
note that the Galactic gas-to-dust ratio assumed above corresponds to a ISM with a higher 
abundance of metals and probably dust grains than that expected to be present in the XUV 
disk of M83, where a larger R g / d might exist (see Galliano et al. 2003, 2005). Moreover, the 
size of the PDRs might be somewhat larger than that of the FUV-bright regions seen in the 
GALEX images, which would lead to a smaller value of Ahi/^4fuv- 

Therefore, we cannot definitively exclude the possibility that a significant fraction of 
the extended HI associated with the XUV disks is a consequence of the photodissociation of 
H 2 in PDRs. 

A more detailed analysis of the relation between gas and recent star formation in XUV 
disks, including a study of the star formation threshold and the star formation law, will be 
the subject of a future paper. 

4.3. Past and future evolution of the XUV disks 

The very low surface brightnesses of these XUV disks at optical wavelengths and, con- 
sequently, the very blue UV-optical colors of their stellar populations led to the suggestion 
that these might be the first generations of stars being formed in these regions (Madore et al. 
2004, Gil de Paz et al. 2005). However, the oxygen abundances derived for the emission-line 
regions of M83 and NGC 4625 (Z~Z /5-Z /lO) indicate that in spite of being relatively 
unevolved chemically these are not the first stars to have formed in the outermost regions of 
these galaxies. Some star formation activity must have taken place in these regions in the 
past in order to enrich the ISM to its present levels. 

There are many ways to produce such an enrichment of the ISM: 

• Low-level continuous star formation: If the star formation history of the XUV disks 
has been similar to that of young spiral disks during their early formation at high 
redshift this level of chemical enrichment could have been reached in a more or less 
continuous way only ~1 Gyr after the first stars formed (see e.g. Boissier & Prantzos 
2000). A scenario where the first stars formed much earlier than one billion years ago 
and have maintained the same level of star formation is unlikely since that would result 
in redder colors than those measured in the XUV disks. In particular, if we correct 
the color profiles of NGC 4625 given in Gil de Paz et al. (2005) for extinction using 
the relation between the (FUV— NUV) color and Apuv found by Boissier et al. (2006) 
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for the disks of spiral galaxies (assuming the Galactic extinction law) the corrected 
(NUV— B) color would range between 0.4 mag (inner XUV disk) and —0.5 mag (outer 
XUV disk). For this range of colors any stellar population with an approximately 
constant star formation rate (and Z=Z Q /5) would necessarily be younger than 1 Gyr 
(e.g. Bruzual & Chariot 2003). On the other hand, the corrected (B — R) color is 
approximately flat with a value of ~0. 6-0. 7 mag. This value is consistent with a 2- 
5 Gyr-old continuously-forming stellar population or with a 500 Myr-old instantaneous 
burst. However, since the photometry errors in (B — R) are large compared with its 
age sensitivity, significantly younger (or older) ages are also possible. 

As noted above, the amount of gas in the XUV disks of these systems is large enough 
to maintain the current level of star formation for several Gyr and up to a Hubble 
time depending on the XUV disk and the region of the XUV disks considered, espe- 
cially once the contribution of a possible molecular component to the total gas mass 
is included. It is therefore conceivable that these galaxies can keep forming stars in 
their outermost regions at the current rate for a few Gyr and that the optical disk will 
steadily grow thanks to the accumulation of the light from the resulting long-lived low- 
and intermediate-mass stars. 

• Episodic star formation: Alternatively, an episodic star formation with brief epochs of 
active star formation followed by long quiescent periods could also have led to this level 
of chemical enrichment. If this were the case, other (now quiescent) galaxies might have 
gone through one or several of these episodes during their lifes. Thilker et al. (2005b) 
have recently proposed that the so-called anti-truncated disks identified by Erwin, 
Beckman, & Pohlen (2005) in early-type spirals could be the evolved counterparts 
of the massive XUV-disk population. The acquisition of deep ground-based multi- 
wavelength surface photometry and single-star color-magnitude diagrams would greatly 
benefit the determination of the star formation history of XUV disks along with its 
possible connection with the anti-truncated disks. Deeper optical spectroscopy of some 
of these regions with the aim of determining more precise O/H and N/ O measurements 
should also provide additional clues on the star formation histories of XUV disks. Some 
of these efforts are currently underway. 

5. Conclusions 



Analysis of the emission-line spectra of a sample of Ha-selected regions in the extended- 
UV disks of M 83 and NGC 4625 allow us to draw the following conclusions regarding the 
properties of these disks: 
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• The line ratios measured in the emission-line regions of the XUV disks of M 83 and 
NGC 4625 show that their ionizing radiation is emitted by young stars and neither due 
to scattering by dust nor to hot evolved stars, such as blue-HB or post-AGB stars. 

• The metal abundances derived, both from the R23 parameter and by comparison of 
multiple line ratios with the predictions of photoionization models, are non primor- 
dial with oxygen abundances between Z /1O and Z /4 in NGC 4625. The oxygen 
abundance in M83 increases from ~Z Q /10 in the outermost regions of the XUV disk 
to close-to-solar (or even super-solar if the M91 calibration is adopted) abundances 
in the optical disk. The same behavior is observed for the radial distribution of the 
dust extinction: it is approximately flat in the case of NGC 4625 and shows a progres- 
sive decrease with increasing galactocentric distance (although with a large scatter) in 
M 83. The comparison of the observed line ratios with the model predictions favors a 
high N/O abundance ratio for the XUV disk of M83. This could be due to the infall 
of pristine gas in regions of secondary nitrogen production or be the natural outcome 
of a top-light IMF. 

• The emission-line luminosities and ratios measured indicate that most of the HII regions 
analyzed are photoionized by single massive stars and only a few are powered by a 
stellar cluster. This fact, along with the low azimuthally- averaged star formation rates 
derived, suggests that stochastic effects both on the upper IMF (Cervino 1998) and 
on the notion of star formation threshold (Boissier et al. 2006) should not be ignored 
when studying the properties of the outer edges of spiral galaxies in general and of 
XUV disks in particular. 

• The amount of gas in these XUV disks is enough to maintain the current level of 
star formation for several Gyr and, in some cases, up to a Hubble time. We cannot 
exclude the possibility that the HI associated with these disks may be the product 
of the dissociation of H 2 in PDRs by the FUV radiation emitted by the stars formed 
during the current and recent episodes of XUV emission. 

• The observed colors and the level of chemical enrichment measured suggest that these 
XUV disks have experienced either (1) a continuous star formation for the last billion 
years or less, or (2) a episodic star formation history with XUV episodes followed by 
quiescent periods. Should the latter scenario be correct, other (perhaps all) spiral 
galaxies might have gone through one or several of these episodes. 

GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in April 2003. 
We gratefully acknowledge NASA's support for construction, operation, and science anal- 
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ysis for the GALEX mission, developed in cooperation with the Centre National d'Etudes 
Spatiales of France and the Korean Ministry of Science and Technology. We thank the 
anonymous referee for his/her constructive comments that have considerably improved the 
content of the paper. AGdP is financed by the MAGPOP EU Marie Curie Research Train- 
ing Network and partially by the Spanish Programa Nacional de Astronomfa y Astrofisica 
under grant AYA2003-01676. We thank Judith Cohen for kindly providing her Ha filter for 
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Table 1. Summary of the imaging observations 
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Table 2. XUV-regions parameters 



XUV region name RA(2000) DEC(2000) f H/3 W,o d 
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f H/3,0 




d 






(h:m:s) 


(d:m:s) 


(ergs -1 cm 


- 2 ) 


(ergs -1 cm 


- 2 ) 


(kpc) 


(1) 




(2) 


(3) 


(4) 




(5) 




(6) 


NGC 4625: 


XUV 02 


12:42:05.136 


+41:16:12.14 


4.10x10" 


16 


1.14x10" 


■15 


6.5 


NGC 4625: 


XUV 03 


12:42:04.025 


+41:16:35.98 


1.73x10' 


■16 


3.71x10" 


■16 


5.9 


NGC 4625: 


XUV 04 


12:42:00.322 


+41:16:39.80 


1.66x10" 


■16 


9.95x10- 


-16 


4.0 


NGC 4625: 


XUV 05 


12:41:59.150 


+41:15:21.11 


2.05x10" 


■16 


3.20x10- 


■16 


4.5 


NGC 4625: 


XUV 06 


12:41:57.029 


+41:15:53.65 


6.75x10- 


■16 


1.57x10" 


•15 


2.7 


NGC 4625: 


XUV 07 


12:41:56.234 


+41:15:30.90 


4.66x10- 


■16 


1.21x10" 


■15 


3.1 


NGC 4625: 


XUV 08 


12:41:54.310 


+41:17:26.82 


1.61x10" 


■16 


7.11x10" 


•16 


2.9 


NGC 4625: 


XUV 09 


12:41:51.538 


+41:15:17.52 


1.02x10" 


■16 


2.78x10- 


■16 


3.2 


NGC 4625: 


XUV 10 


12:41:48.677 


+41:17:46.81 


1.01x10" 


15 


1.01x10" 


-15 


4.4 


NGC 4625: 


XUV 11 


12:41:45.151 


+41:16:00.29 


2.60x10" 


■16 


2.35x10- 


-15 


4.3 


NGC 4625: 


XUV 12 


12:41:42.977 


+41:16:37.95 


5.33x10- 


■16 


1.10x10" 


•15 


3.7 



Note. — (1) Region identification, (2) Right Ascension, (3) Declination, (4) Observed H/3 
flux inside the corresponding slitlet. (5) Extinction-corrected H/3 flux inside the slitlet; color 
excesses are given in Table 3. (6) Galactocentric distance in kpc. '''Region belongs to the 
optical disk of M83. 



Table 3. XUV-regions line ratios 



XUV region name E(B - V) [OIII]/[OII] [OIII]/H/3 [NII]/Hq [SII]/Hq 6717/6731 R23 log(EW HQ ) log(EW[ OIII ] ) Z 1/h M91 Z 1/h pl 
(!) ( 2 ) (3) ( 4 ) (5) (6) (7) (8) (9) (10) (11) (12) 



M83 (NGC 5236) 



NGC 


5236 


XUV 


01 





00 







38 


0.49 


-0 


89 


-1.11 


1 


79 





76 


2 


71 


2 


28 


7.861/8 


680 


7.718/8.458 


NGC 


5236 


XUV 


02 





23 




>0 


76 


0.86 


-1 


24 


<-0.73 






<1 


05 


>1 


80 


>1 


91 








NGC 


5236 


XUV 


03 





29 


> 


-0 


20 


0.14 


-0 


65 


<-0.47 






<0 


68 


2 


89 


2 


12 








NGC 


5236 


XUV 


04 





00 








<-0.21 


-0 


66 


<-0.11 






<0 


30 


2 


23 


<2 


00 








NGC 


5236 


XUV 


05 





53 


< 


-0 


42 


<-0.06 


-0 


37 


-0.50 


1 


04 


<0 


63 


2 


51 


<3 


10 








NGC 


5236 


XUV 


06 





00 


< 


-0 


52 


<0.01 


-0 


51 


-0.65 





88 


<0 


76 


2 


96 


<1 


53 








NGC 


5236 


XUV 


07 





44 




-0 


73 


-0.07 


-0 


60 


-0.98 


2 


40 





85 


2 


56 


1 


76 


8.391/8 


427 


7.765/7.744 


NGC 


5236 


XUV 


08 





07 


< 


-0 


78 


<0.12 






<-0.19 






<1 


09 


2 


75 


<1 


64 








NGC 


5236 


XUV 


09 





12 







06 


0.39 


-0 


89 


-0.88 


2 


29 





79 


3 


38 


2 


67 


8.015/8 


617 


7.891/8.316 


NGC 


5236 


XUV 


10 





10 




-0 


28 


0.03 


-0 


81 


-0.76 


1 


11 





62 


2 


99 


2 


10 


7.862/8 


761 


7.506/8.341 


NGC 


5236 


XUV 


11 





00 







17 


0.48 


-1 


02 


-1.11 


2 


01 





83 


2 


52 


2 


42 


8.049/8 


584 


7.952/8.311 


NGC 


5236 


XUV 


12 





49 




-0 


02 


0.16 


-0 


52 


-0.84 


1 


07 





59 


2 


72 


2 


34 


7.730/8 


798 


7.550/8.492 


NGC 


5236 


XUV 


13 





19 







04 


0.32 


-0 


69 


-1.01 


1 


21 





72 


3 


39 


3 


29 


7.904/8 


689 


7.751/8.391 


NGC 


5236 


XUV 


14 





34 




-0 


29 


-0.13 


-0 


46 


-0.74 





87 





47 


2 


63 


1 


63 


7.660/8 


875 


7.306/8.468 


NGC 


5236 


XUV 


15 





00 


> 


-0 


56 


-0.66 


-0 


39 


-0.60 





77 


<0 


14 


2 


88 


1 


90 








NGC 


5236 


XUV 


16 





28 


< 


-0 


78 


<-0.46 


-0 


34 


-0.72 


1 


60 


<0 


51 


>1 


91 


<0 


81 








NGC 


5236 


XUV 


17 





13 




-0 


58 


0.05 


-0 


85 


-0.65 





83 





86 


3 


03 


1 


93 


8.348/8 


442 


7.877/7.824 


NGC 


5236 


XUV 


18 





18 








0.13 


-0 


48 


-0.93 


1 


06 






2 


85 


3 


10 








NGC 


5236 


XUV 


19 





18 








<0.03 


-0 


43 


-0.31 





97 


<0 


74 


>1 


41 


<1 


84 








NGC 


5236 


XUV 


20 





41 




-0 


49 


-0.54 


-0 


38 


<-1.82 









20 


2 


76 


1 


23 


... /9 


007 


• • • /8.543 


NGC 


5236 


XUV 


21 





39 




-0 


29 


-0.59 


-0 


38 


-0.93 


1 


31 





01 


3 


05 


1 


78 


... /9 


061 


• • • /8.696 


NGC 


5236 


XUV 


22 





27 







02 


-0.27 


-0 


38 


-0.78 


1 


30 





14 


>2 


67 


3 


54 


... /9 


026 


• • • /8.759 


NGC 4625 


NGC 


4625 


XUV 


01 





31 







51 


0.59 


-1 


26 


<-0.59 









83 


2 


85 


2 


54 


7.941/ 




7.787/ ■ ■ • 


NGC 


4625 


XUV 


02 





31 




-0 


43 


0.30 


-1 


05 


-0.94 


1 


90 





99 


2 


76 


2 


24 


8.531/ 




8.486/ • • • 


NGC 


4625 


XUV 


03 





23 








<-0.23 






-0.51 


1 


63 


<0 


69 


2 


82 


<2 


28 








NGC 


4625 


XUV 


04 





54 


< 


-0 


99 


<-0.23 


-0 


91 


-0.63 


1 


19 


<0 


92 


2 


77 


<1 


65 








NGC 


4625 


XUV 


05 





13 




-0 


25 


0.22 


-0 


81 


-0.79 





90 





79 


2 


37 


2 


00 


8.123/ 




7.868/ • ■ • 


NGC 


4625 


XUV 


06 





25 




-0 


65 


-0.07 


-0 


81 


-0.66 





68 





79 


2 


57 


1 


68 


8.265/ 




7.635/ ■ ■ • 


NGC 


4625 


XUV 


07 





29 




-0 


09 


0.44 


-1 


09 


-1.04 


1 


20 





92 


>2 


32 


>2 


07 


8.296/ 




8.269/ ■ ■ • 


NGC 


4625 


XUV 


08 





45 








0.52 


-0 


99 


<-0.69 










3 


22 


2 


63 








NGC 


4625 


XUV 


09 





30 







25 


0.96 


-0 


85 


<-0.44 






1 


28 


3 


02 


3 


10 


8.400/ 




8.400/ ■ • • 


NGC 


4625 


XUV 


10 





00 








0.68 


-1 


29 


-0.66 


1 


36 






3 


12 


2 


95 








NGC 


4625 


XUV 


11 





66 




>0 


26 


0.92 


-1 


18 


<-0.68 






<1 


24 


2 


61 


>2 


04 








NGC 


4625 


XUV 


12 





22 




-0 


55 


0.07 


-1 


10 


-0.69 


1 


73 





85 


3 


07 


2 


32 


8.333/ 




7.891/ ••• 



CO 



Note. — (1) Region identification, (2) color excess derived from the Hq/H/3 Balmcr-linc decrement; it includes the effects of both the Galactic and the 
internal dust extinction, (3) logarithm of the extinction-corrected [OIII] A5007A/[OII] AA3726.3729AA line ratio, (4) logarithm of the [OIII] A5007A/H/3 line 
ratio, (5) logarithm of the [Nil] A6584A/Ha line ratio, (6) logarithm of the [SII] AA6717,673lAA/Ha line ratio, (7) [SII] A6717A/ [SII] A6731 A line ratio, (8) 
R23 parameter, defined as log{ ([Oil] AA3726,3729AA+[OIII] A5007A+[OTII] A4959A)/H/3} , (9) logarithm of the equivalent width (in angstroms) of Ha, (10) 
logarithm of the equivalent width of [OIII] A5007A, (11) nebular oxygen abundance, 12 + log(0/H), for the low/high-mctallicity branch of the R23 calibration 
using the method of McGaugh (1991; M91), (12) oxygen abundances for the low/high- metallicity branch of the R23 calibration using the method of Pilyugin 
& Thuan (2005; PT05). We have adopted 12+log(O/H) =8.69 throughout the text (Allcndc Pricto ct al. 2001). 



Table 4. Best-fitting single-star photoionization models (solar N/O ratio) 



XUV region name \ 2 



(1) 



(2) 



fH/3,modol 

(ergs" 1 cm~ 2 ) 
(3) 



Age* 
(Myr) 
(4) 



T e ff,* 

(K) 
(5) 



log(g)* 
(6) 



Mass, 
(M Q ) 
(7) 



^Jgas 

(Z ) 
(8) 



log(n gas ) 
(cm -3 ) 
(9) 



l0g(Rin 

(cm) 
(10) 



M83 (NGC 5236) 



NGC 


5236: 


XUV 


01 





.355 


1 


.207E- 


i rl.207B- 
D.lbuii — 


■15 
■16 


3 


A.3.0 

O.U 




3. 


eo3.68 

O.Do 


40. 


n.40.o 

40 


a oaO.60 
u - zu 0.20 


2. 


30 

2 


17. 


cl7.5 

- O Q 


NGC 


5236: 


XUV 


07 


1 


.975 


3 


.196E- 


1fi 3.203B- 
1U 3.052_E- 


16 
16 


3 


a3.0 
• u 3.0 


35636I1IH 


3. 


n fi 3.96 
yD 3.96 


25. 


A25.0 
u 25.0 


a oaO.30 
U.ZUq.20 


2. 


A2.0 
u 2.0 


16. 


A17.5 
• u 16.0 


NGC 


5236: 


XUV 


09 





.383 


4 


.321E- 


-, pAAWE- 
iD 4.117_E- 


15 
16 


1 


n 30 


0770740576 
01 '"'36625 


4. 


it;4.15 
10 3.68 


25. 


A.60.0 
u 25.0 


90°- 3 
U.ZUq 20 


2. 


A2.0 
U 0.0 


17. 


A17.5 
• u 16.0 


NGC 


5236: 


XUV 


10 





.413 


1 


.379E- 


-, rl.379-E- 

' 10 3.196.E- 


15 
16 


3 


n3.0 
U 3.0 


ODDZ0 356 3 6 


3 


R o3.96 
D8 3.68 


40. 


a40.0 
u 25.0 





in - 20 

1U 0.10 


0. 


n 20 


17. 


cl7.5 
°16.0 


NGC 


5236: 


XUV 


11 





.828 


1 


.207E- 


-, rl.207B- 
lo 1.201.E- 


15 
15 


3 


n 30 


36625 3 gg25 


3 


eo3.68 
D8 3.68 


40. 


A40.0 
■ u 40.0 





oa0.20 
• zu 0.20 


2. 


a2.0 
u 2.0 


17. 


C17.5 
°16.0 


NGC 


5236: 


XUV 


12 





.427 


2 


.592E- 


1fi 2.733B- 
" 1U 2.385.E- 


16 
16 


3 


o 30 


^fi^tf; 35636 
ooooo 35636 


3 


n fi 3.96 
yD 3.96 


25. 


A25.0 
u 25.0 





fin - 80 

OU 0.40 


3. 


A3.0 
U 3.0 


17. 


C17.5 
°17.5 


NGC 


5236: 


XUV 


13 





.207 


2 


.677E- 


1f; 2.901B- 
1D 2.326.E- 


16 
16 


3 


n 30 

U 3.0 


35636 3 5g3g 


3. 


n R 3.96 
yD 3.96 


25. 


A25.0 
u 25.0 





oa0.50 
• du 0.20 


3 


p3.0 
U 3.0 


17. 


A17.5 
• u 16.0 


NGC 


5236: 


XUV 


14 





.281 


1 


.654E- 


1fi 1.654B- 
10 1.385.E- 


16 
16 


1 


• u 1.0 


34984IH81 


4. 


1 «4.18 
io 4.18 


20. 


A20.0 
• u 20.0 





OAl.OO 

•oU .80 


2. 


A2.0 
u 2.0 


16. 


A16.0 
u 16.0 


NGC 


5236: 


XUV 


17 


1 


.848 


1 


.351E- 


-, rl.351.E- 
10 1.351,E- 


15 
15 


3 


n 30 


'iRR'? C36625 


3 


eo3.68 
°°3.68 


40. 


A40.0 
■ u 40.0 





9A0.20 


0. 


aO.O 
U 0.0 


16. 


C17.5 
°16.0 


NGC 


5236: 


XUV 


20 





.419 


1 


.044E- 


i fi 1.044_E- 
iD 6.723B- 


16 
17 


1 


n 30 


34984 3 f 7 lf 


4. 


1 «4.18 
• io 4.05 


20. 


A20.0 
• u 20.0 


1 


A A 1.50 
UU 1.00 


3. 


A3.0 
U 3.0 


17 


C17.5 
°16.0 


NGC 


5236: 


XUV 


21 





.741 


9 


.366E- 


i 7 9.366£- 
1 '5.452_E- 


17 
17 


1 


o 30 


34984 3 f 7 lf 


4. 


1 «4.18 
■ io 4.05 


20. 


A20.0 
• u 20.0 


1 


CA2.00 
OU 1.50 


3. 


A3.0 
U 3.0 


17. 


C17.5 
°16.0 


NGC 


5236: 


XUV 


22 


1 


.486 


6 


.525E- 


-, 7 7.187.E- 
1 '6.362E- 


17 
17 


1 


n i.o 

• u 1.0 


34984I1H1 


4. 


1 «4.18 
404.18 


20. 


A20.0 
• u 20.0 


1. 


ca1.50 
• ou 1.50 


3 


a3.0 
u 3.0 


16. 


C17.0 
•°16.0 



NGC 4625 



NGC 4625: XUV 01 0.702 



1.126E-15| : |«i = i| 



NGC 4625: XUV 02 4.260 4.640E-16 1 ; I 3 ?I: 1 6 



l-0?;8 
i-o?:8 



45087fH|I 
42600| 2 i° 



3.99 3 ;i 
4.08l°| 



60.0| 20 ° 

40.0^:8 



o-io8;?8 
o-io8:i8 



1 n 30 
n o 10 



lo - u 16.0 
lo -°16.0 



Table 4 — Continued 



XUV region name 
(1) 


2 

X 
(2) 


fH/3, model 

(ergs" 1 cm" 2 ) 
(3) 




* 

Age* 
(Myr) 
(4) 


T e ff,* 

(K) 
(5) 


log(g)* 
(6) 


Mass* 
(M ) 
(7) 


^gas 

(Zo) 
(8) 


log(n gas ) 
(cm" 3 ) 
(9) 


1 /T~> \ 

l0g(Ri„) 

(cm) 
(10) 


NGC 4625: XUV 05 


0.609 


2.914E-16 2 Sl 


16 
16 


3-Of;8 


^tfifiC^ 37707 


3-683:68 


40.oi:° 


9f)0-30 


1 n 10 

LM 1.0 


1fi o 17 - 5 

1D - U 16.0 


NGC 4625: XUV 06 


1.478 


3.031E-16?;°^|l 


16 
16 


3-0?;8 


^fifiC^ 37707 


3-68|;68 


4o.oi:« 


9n°- 20 

u - zu 0.10 


n n°° 
u - u o.o 


i '- u 16.0 


NGC 4625: XUV 07 


2.504 


2-695E-16^f: 


15 
16 


3-Of;8 


36625!<3 (325 


3-683:68 


40.01°:" 


90 - 20 


9 n2.0 


1 a rl7.5 
1D -°16.0 


NGC 4625: XUV 09 


5.162 


3.218E-15^fl 


15 
15 


1 O 10 


49565|ii 


093-92 
°- az 3.92 


120.0118:8 


o.eo":!" 


1 o 10 
1 - u o.o 


io - u 16.0 


NGC 4625: XUV 12 


3.405 


1.056E-16?Sl 


16 
16 


1 n 30 


37707 3 I 6 7 i 


4 - i0 3.68 


25.oi:° 


0-108:18 




1 7 n 17 - 5 

1 ' - u 16.0 



Note. — (1) Region identification, (2) best-fitting reduced x 2 , (3) H/3 flux predicted by the photoionization models, (4) age of 
the star (Myr), (5) effective temperature (in K), (6) logarithm of the gravity, (7) ZAMS mass (M ) , (8) metal abundance of the 
gas and dust grains (in Z ), (9) density of the gas (in atoms cm" 3 ), (10) inner radius of the spherical gas distribution (in cm). 
The 68.3% confidence intervals for each parameter considered separately are also given. 



Table 5. Best-fitting single-star photoionization models (O/H-scaled N/O ratio) 



XUV region name x 2 W.modei Age* T effj * log(g)* Mass* Z gas log(n gas ) log(R in ) 

(ergs" 1 cm- 2 ) (Myr) (K) (M Q ) (Z Q ) (cm- 3 ) (cm) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



M83 (NGC 5236) 



NGC 


5236: 


XUV 


01 


3 


.022 


7. 


.773E-16^fl 


-15 

- IO 


1 


nl.0 
,u 1.0 


42600|f^f 


4. 


08||g2 


40.0i|° ° 


2 


■001:8!] 


2 


n3.0 

■Ul_Q 


16. 


.0 


NGC 


5236: 


XUV 


07 


5 


.030 


1 


.537E-16};^|: 


16 
16 


1 


nl.O 
u 1.0 


34984|t||| 


4. 


io4.18 
i8 4.18 


20.ol8;8 


1 


nnl.00 
uu 1.00 


2 


n2.0 
• u 2.0 


16. 


.0 


NGC 


5236: 


XUV 


09 


2 


.892 


1 


.424E-14^ = 


14 
■16 


1 


n 30 

u 1.0 


49565I11H 


3. 


003.96 
yz 3.92 


120.0^0° 


1 


cnl.50 
ou 1.00 





n3.0 
• u 0.0 


17. 


.0 


NGC 


5236: 


XUV 


10 





.905 


1 


■685E-16^f = 


16 
-16 


1 


nl.O 
• U 1.0 


34984 3 !i| 


4. 


iS 4.18 


on r,20.0 
ZU - U 20.0 





on0.80 
•° U 0.80 


2 


• u 2.0 


16. 


.0 


NGC 


5236: 


XUV 


11 


2 


.589 


1 


,716E-16 2 ;^|: 


-16 
-16 


3 


n 30 

U 3.0 


3563635535 


3. 


Ofi3.96 
yD 3.96 


oc A25.0 
ZO - U 25.0 


1 


nnl.00 
UU 0.80 


3 


n 30 

• U 3.0 


16. 


.5 


NGC 


5236: 


XUV 


12 


6 


.731 


1 


.685E-16S: 


16 
■16 


1 


nl.O 
U 1.0 


34984il|H 


4. 


io4.18 
i8 4.18 


on n 20 

ZU - U 20.0 





onl.00 
•° U 0.80 


2 


o 20 


16. 


.0 


NGC 


5236: 


XUV 


13 


5 


.489 


1 


.461E-16l;^ = 


14 
16 


3 


n 30 

u 1.0 


35636|11|1 


3. 


96l;l| 


oc nl20.0 
ZO - U 25.0 


1 


cnl.50 
■ ou 1.00 


3 


• u o.o 


17. 


.0 


NGC 


5236: 


XUV 


14 


4 


.918 


6 


.495E-17j;^I 


16 
17 


1 


nl.O 
U 1.0 


34984I1H1 


4. 


io4.18 
i8 4.18 


zu - u 20.0 


1 


cnl.50 
ou 1.00 


3 


• U 2.0 


16. 


.5 


NGC 


5236: 


XUV 


17 


1 


.338 


1 


■685E-16};^|: 


16 
-16 


1 


nl.O 
U 1.0 


34984I1H1 


4. 


104.18 
i8 4.18 


on n20.0 
ZU - U 20.0 





on0.80 
•° U 0.80 


2 


n 20 

• U 2.0 


16. 


.0 


NGC 


5236: 


XUV 


20 


5 


.418 


6 


,160E-17^23| : 


-17 
-17 


1 


n 30 

u 1.0 


34984|^84 


4. 


io4.18 
i8 4.05 


on n20.0 
ZU - U 20.0 


2 


nn 200 

• uu 1.50 


3 


n 30 

• U 3.0 


17. 


.0 


NGC 


5236: 


XUV 


21 


3 


.836 


4 


.705E-17S: 


■17 
17 


3 


n3.0 
u 1.0 


007C1 34984 
°°' 01 33751 


4. 


uo 4.05 


on n 20 

ZU - U 20.0 


2 


oo 200 

uu 2.00 


3 


O 30 


16. 


.5 


NGC 


5236: 


XUV 


22 


4 


.836 


5 


.465E-17«;^|: 


17 
17 


1 


nl.O 
U 1.0 


34984I1H1 


4. 


io4.18 
i8 4.18 


ZU - U 20.0 


2 


nn 200 

UU 2.00 


3 


n3.0 
• U 3.0 


16. 


.5 



17.5 
16.0 

16.0 
16.0 



,16.0 
16.0 



16.0 
16.0 



,16.0 
16.0 



17.5 
16.0 



NGC 4625 

NGC 4625: XUV 01 0.117 7.990E-16 2 .;*™! l.O^g 45087^1? 3-99|;f§ 60. O^ l.OO^gg 1.0j;8 17.5 
NGC 4625: XUV 02 4.120 5.735E-17^^|l^ 5.0f;g 31698|f|g| 3.65|;^ 25.0 2 8 : 8 0.50g : fg l.Of g 16.0 



Table 5 — Continued 



XUV region name 
(1) 


9 

X 
(2) 


fH/3, model 

(ergs -1 cm -2 ) 
(3) 




Age* 
(Myr) 
(4) 


T e ff,* 

(K) 
(5) 


log(g)* 
(6) 


Mass* 
(M ) 
(?) 


Zgas 

(Z ) 
(8) 


log(n gas ) 
(cm -3 ) 
(9) 


1 /T~l \ 

l0g(Ri n ) 

(cm) 
(10) 


NGC 4625: 


XUV 05 


2.504 


3-781E-17^f : 


17 
17 


i nl.O 
L - u 1.0 


349S4it||t 


4.18HI 


on n20.0 
zu - u 20.0 


0-808:18 


9 2 
z - u l.O 


i6.oll : 8 


NGC 4625: 


XUV 06 


1.256 


3.781E-17^ifl 


17 
17 


i nl.O 
L - u 1.0 


34984it||t 


4.18HI 


90 20 

zu - u 20.0 


0-808:18 


9 2 

z - u 2.0 


ll:) . u 16.0 


NGC 4625: 


XUV 07 


2.333 


4.927E-17^fl 


17 
17 


1 n 10 

L - u 1.0 


34984|tl|l 


4.18HI 


90 20 

zu - u 20.0 


0.508:18 


9 2 

z - u 2.0 


iD - u 16.0 


NGC 4625: 


XUV 09 


0.280 


1.839E-16 i; ^| : 


16 
16 


i nl.O 
L - u 1.0 


34984it||t 


4.181H 


90 20 

zu - u 20.0 


0.80818 


O 


1D - 16.5 


NGC 4625: 


XUV 12 


1.366 


3.781E-17|^fl 


17 
17 


i nl.O 


34984it||t 


4.18HI 


on n20.0 
ZU - U 20.0 


0.808I8 


9 n 2 

Z - U 2.0 


ifi n 160 

lo - u 16.0 



Note. — (1) Region identification, (2) best-fitting reduced x 2 , (3) H/3 flux predicted by the photoionization models, (4) age 
of the star (Myr), (5) effective temperature (in K), (6) logarithm of the gravity, (7) ZAMS mass (M ) , (8) metal abundance of 
the gas and dust grains (in Z ), (9) density of the gas (in atoms cm -3 ), (10) inner radius of the spherical gas distribution (in 
cm). The 68.3% confidence intervals for each parameter considered separately are also given. 
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Table 6. Best-fitting starburst photoionization models (solar N/O ratio) 



XUV region name x 2 iH/3,modei Age starburst Z gas n gas log(R in ) 

(ergs _1 cm~ 2 ) (Myr) (Z ) (cm~ 3 ) (cm) 

(1) (2) (3) (4) (5) (6) (7) 



M83 (NGC 5236) 



NGC 5236: 


XUV 


01 


1.845 


7.972E- 


ql.022£- 
' L °7.972E- 


12 
13 


r. n5.0 
°- u 5.0 


n 2n - 20 


1 rl.5 

i -°1.0 


iD -°16.0 


NGC 5236: 


XUV 


07 


2.487 


6.075E- 


i o7.171£- 
io 6.075E- 


13 
13 


6 n 60 


n 2n - 30 

u - zu 0.20 


1 O 10 


10.0 16 Q 


NGC 5236: 


XUV 


09 


2.471 


1.022E- 


i 9 1.238£- 
iz 1.022£- 


12 
12 


°- U 5.0 


20 20 

u - zu 0.20 


1 L0 


16 17 - 5 


NGC 5236: 


XUV 


10 


2.414 


1.238E- 


i 9 1.238£- 
±Z 2.515E- 


12 
13 


r n5.0 
°- u 5.0 


20 60 

u - zu 0.20 




If! c;17.5 


NGC 5236: 


XUV 


11 


1.776 


8.425E- 


i a7.972E- 


13 
-14 


2 5 - 


i crnl.50 
L - ou 0.20 


2.0?:^ 


16 n 17 - 5 


NGC 5236: 


XUV 


12 


3.387 


1.238E- 


-, 9 1.238£- 
iz 3.695£- 


12 
13 


°- u 5.0 


n 2n - 40 

u - zu 0.20 


u -°0.5 


17 n 17 - 5 

1 ' - u 16.0 


NGC 5236: 


XUV 


13 


1.372 


1.238E- 


-, 9 1.238£- 
iz 1.238£- 


12 
12 


r. n5.0 
°- u 5.0 


n 2n - 20 

u - zu 0.20 


u -°0.5 


16 n 17 - 5 

iu - u 16.0 


NGC 5236: 


XUV 


14 


2.488 


6.754E- 


i o7.899.E- 
i0 6.754E- 


13 
13 


r. n5.0 
°- u 5.0 


n 4n - 40 

u -^ u 0.30 


o s - 5 

u - d 0.5 


1 ' -°16.0 


NGC 5236: 


XUV 


17 


4.780 


1.238E- 


i 9 1.238£- 
iz 1.238£- 


12 
12 


r q5.0 
°- u 5.0 


20 - 20 

u - zu 0.20 


u -°0.5 


17 17 - 5 

1 ' - u 16.0 


NGC 5236: 


XUV 


20 


2.350 


1.565E- 


i ol.852.E- 
io 7.287.E- 


13 
15 


6 6 

u - u 6.0 


60 L0 ° 


1 c3.5 
L -°1.5 


1fi ^i 17 - 5 


NGC 5236: 


XUV 


21 


4.106 


1.572E- 


-i ol.572.E- 
io 1.565E- 


13 
13 


6 n 60 

°- u 6.0 


n 6n - 60 

u - uu 0.60 


1 "S 1 - 5 
L -°1.5 


1 ' - d 16.0 


NGC 5236: 


XUV 


22 


4.888 


5.207E- 


i o5.896E- 
io 5.207E- 


13 
13 


r. n5.0 
°- u 5.0 


n 6n - 60 


u -°0.5 


17 n 17 - 5 

1 ' - u 16.0 



NGC 4625 



NGC 4625: 


XUV 01 


0.519 


1.020E-14™: 


-13 

-15 


2 5 - 


i ^n 2 00 


o c3.5 
z - d 0.5 


17.5 


NGC 4625: 


XUV 02 


5.265 


1.890E-14 2 ; 8 7 ^: 


-13 
-14 


2 5 - 


i - du 0.20 


2 n 2 

z - u 0.5 


16.0 


NGC 4625: 


XUV 05 


2.255 


2.777E-13 2 ;^I 


-13 
-13 


c: n5.0 
°- u 5.0 


u - zu 0.20 


u -°0.5 


16.0 


NGC 4625: 


XUV 06 


4.145 


1.405E-13}; 77 , 2 !: 


-13 

-13 


r. n5.0 
°- u 5.0 


n 3n - 40 


1 o LO 


16.0 


NGC 4625: 


XUV 07 


2.453 


1.890E-14JS: 


-14 
-14 


2 n 20 


i - ou 1.50 


2 n 2 - 


16.0 



17.5 
16.0 

17.5 
16.0 

17.5 
16.0 

17.5 
16.0 

17.5 
16.0 
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Table 6 — Continued 



XUV region name 
(1) 


x 2 

(2) 


lH/3,model 

(ergs -1 cm -2 ) 
(3) 




Age s tarburst 

(Myr) 
(4) 


Zgas 

(z ) 

(5) 


ftgas 

(cm -3 ) 
(6) 


l0g(Rin) 

(cm) 
(7) 


NGC 4625: XUV 09 


8.609 


2.758E-15^|: 


-15 
-15 


2.ol;8 


2.ool;88 


3.o|8 


16.0H8 


NGC 4625: XUV 12 


6.202 


3.644E-14*;™: 


-13 
-14 


2.0|8 


i-505:fg 


i-55:I 


i6.5^ : 8 



Note. - (1) Region identification, (2) best-fitting reduced \ 2 ■, (3) H/3 flux predicted by 
the photoionization models, (4) age of the starburst (Myr), (5) metal abundance of the gas 
and dust grains (in Z ), (6) density of the gas (in cm -3 ), (7) inner radius of the spherical gas 
distribution (in cm). A stellar mass of 1O 3 M was adopted for the ionizing starburst. The 
68.3% confidence intervals for each parameter considered separately are also given (see text 
for more details). 
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Table 7. Best-fitting starburst photoionization models (O/H-scaled N/O ratio) 



XUV region name x 2 fH/3,modei Age s tarburst Z gas Ugas log(Rj n ) 

(ergs _1 cnT 2 ) (Myr) (Z Q ) (cm -3 ) (cm) 

(1) (2) (3) (4) (5) (6) (7) 



M83 (NGC 5236) 



NGC 5236: 


XUV 


01 


3.148 


5.316E- 


-i a 5A90E- 
' 1 ^5.316.E- 


•14 
•14 


2 2 - 


2 


OO 200 


2 2 - 


1 6 17 - 5 


NGC 5236: 


XUV 


07 


15.092 


7.254E- 


i P-7.254E- 
io 7.254£- 


■15 
•15 


6 6 


1 


on 100 

• uu 1.00 




16 160 


NGC 5236: 


XUV 


09 


4.676 


5.316E- 


i A5A90E- 
L ^b.3lGE- 


•14 
-14 


2 20 


2 


QQ2.00 


2 2 - 


16 17 - 5 


NGC 5236: 


XUV 


10 


8.301 


1.066E- 


-i ql.073.E- 
lo 1.063E- 


13 
•13 


2 2 - 


2 


QQ2.00 


i -°1.5 


17 17 - 5 

1 ' - u 16.0 


NGC 5236: 


XUV 


11 


1.969 


5.316E- 


-i A5A90E- 
L ^b.3lGE- 


•14 
■14 


2 2 - 


2 


oo 2 - 00 

uu 2.00 


2 2 


iu - u 16.0 


NGC 5236: 


XUV 


12 


14.941 


5.316E- 


i A 5A90E- 
Lq: 5.316E- 


-14 
■14 


2 2 - 


2 


oo 2 - 00 


2 2 - 


16 17 - 5 


NGC 5236: 


XUV 


13 


7.870 


5.316E- 


-i A 5A90E- 
Lq: 5.316E- 


-14 
-14 


2 2 - 


2 


oo 2 - 00 

uu 2.00 


2.02:8 


1 6 17 - 5 


NGC 5236: 


XUV 


14 


15.620 


3.038E- 


-i p-3.038£- 
id 3.038£- 


-15 
•15 


c n5.0 
d - u 5.0 


2 


oo 2 - 00 


3.5l:| 




NGC 5236: 


XUV 


17 


11.446 


1.613E- 


1 ol.628.E- 
io 1.613E- 


-13 
13 


2 2 - 


1 


cnl.50 
ou 1.50 




1 6 17 - 5 


NGC 5236: 


XUV 


20 


3.247 


3.038E- 


-i P-3.038E- 
±J 3.038£- 


-15 
-15 


c n5.0 
°- U 5.0 


2 


QQ2.00 


q cr3.5 

.°3.5 


17 17 - 
i i .u 17 Q 


NGC 5236: 


XUV 


21 


7.273 


3.038E- 


-i P-3.038E- 
± < J 2.441E- 


-15 
-15 


c n5.0 
°- u 5.0 


2 


oo 2 - 00 

• uu 2.00 


3-5|:| 


1 7 17 - 

1 ' - u 16.0 


NGC 5236: 


XUV 


22 


14.128 


3.038E- 


1 p-3.038E- 
id 3.038E- 


-15 
•15 


c n5.0 
°- u 5.0 


2 


oo 2 - 00 

• uu 2.00 


3-5|J 


1 7 n 17 - 

1 ' - u 17.0 



NGC 4625 



NGC 4625: 


XUV 01 


0.640 


6.347E-15^ 7 |: 


•15 
•15 


2 2 - 


2 OO 2 00 

uu 2.00 


2.51J 


17 S 17 - 5 

1 ' - o 16.0 


NGC 4625: 


XUV 02 


6.409 


1.193E-14} : 23 , 2 f: 


-14 
-14 


2 2 


2 OO 2 00 

z - uu 2.00 


2.0|8 


iu . u 16.0 


NGC 4625: 


XUV 05 


8.382 


1.193E-14}; 232 !: 


-14 
14 


2 2 - 

z - u 2.0 


2.00 2 ;™ 


2 2 

z " u 2.0 


16 17 - 5 


NGC 4625: 


XUV 06 


12.131 


2.385E-14K£fl 


■14 
■14 


2 2 


2.00 2 ;™ 


i -°1.0 


16 17 - 5 


NGC 4625: 


XUV 07 


3.150 


1.193E-14}; 232 !: 


-14 
-14 


2 2 


2.oo 2 : °° 


2 2 - 

z - u 2.0 


16 17 - 5 
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Table 7— Continued 



XUV region name 
(1) 


x 2 

(2) 


lH/3,model 

(ergs -1 cm -2 ) 
(3) 




Age s tarburst 

(Myr) 
(4) 


Zgas 

(z ) 

(5) 


Ugas 

(cm -3 ) 
(6) 


l0g(Rin) 

(cm) 
(7) 


NGC 4625: XUV 09 


7.502 


2.758E-15f^!: 


-15 
-15 


2.02:8 


2.00I88 


3.o|8 


16-0118 


NGC 4625: XUV 12 


6.645 


3.619E-14|^|- 


-14 
-14 


2.0|8 


1.50H8 


i-5l:I 


16-0&8 



Note. - (1) Region identification, (2) best-fitting reduced \ 2 ■, (3) H/3 flux predicted by 
the photoionization models, (4) age of the starburst (Myr), (5) metal abundance of the gas 
and dust grains (in Z ), (6) density of the gas (in cm -3 ), (7) inner radius of the spherical gas 
distribution (in cm). A stellar mass of 1O 3 M was adopted for the ionizing starburst. The 
68.3% confidence intervals for each parameter considered separately are also given (see text 
for more details). 



Table 8. Observed and best-fitting line ratios 



XUV region name Observed line ratios Single-star, solar-N/O best fit Single-star, O/H-scalcd-N/O best fit Starburst, solar-N/O best fit Starburst, O/H-scalcd-N/O best fit 

OII/H/3 OIII/H/3 NII/H/3 SII/H/3 OII/H/3 OIII/H/3 NII/H/3 SII/H/3 OII/H/3 OIII/H/3 NII/H/3 SII/H/3 OII/H/3 OIII/H/3 NII/H/3 SII/H/3 OII/H/3 OIII/H/3 NII/H/3 SII/H/3 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (IS) (19) (20) (21) 



M83 (NGC 5236) 



NGC 


5236: 


XUV 


01 


+ 


110 


+0 


490 


-0 


434 


-0 


654 


+ 


116 


+0 


410 


-0 


484 


-0 


509 


-0 


172 


+0 


476 


-0 


623 


-0 


209 


+ 


122 


+0 


362 


-0 


521 


-1 


100 


+ 


116 


+0 


397 


-0 


738 


-0 


968 


NGC 


5236: 


XUV 


07 


+0 


660 


-0 


070 


-0 


144 


-0 


524 


+0 


291 


+0 


019 


-0 


215 


-0 


377 


+0 


259 


-0 


182 


-0 


449 


-0 


216 


+0 


327 


-0 


179 


-0 


205 


-0 


895 


+0 


339 


-0 


550 


-0 


612 


-1 


182 


NGC 


5236: 


XUV 


09 


+0 


330 


+0 


390 


-0 


434 


-0 


424 


+0 


146 


+0 


393 


-0 


422 


-0 


407 


+0 


011 


+0 


285 


-0 


606 


-0 


077 


+0 


198 


+0 


263 


-0 


414 


-0 


970 


+0 


116 


+0 


397 


-0 


738 


-0 


968 


NGC 


5236: 


XUV 


10 


+0 


310 


+0 


030 


-0 


354 


-0 


304 


+0 


216 


-0 


067 


-0 


408 


-0 


338 


+0 


364 


-0 


028 


-0 


520 


-0 


182 


+0 


240 


+0 


176 


-0 


348 


-0 


844 


-0 


041 


-0 


027 


-0 


739 


-1 


013 


NGC 


5236: 


XUV 


11 


+0 


310 


+0 


480 


-0 


564 


-0 


654 


+0 


116 


+0 


410 


-0 


484 


-0 


509 


+0 


043 


+0 


277 


-0 


709 


-0 


456 


+0 


204 


+0 


387 


-0 


751 


-0 


954 


+0 


116 


+0 


397 


-0 


738 


-0 


968 


NGC 


5236: 


XUV 


12 


+0 


180 


+0 


160 


-0 


064 


-0 


384 


+0 


121 


+0 


232 


-0 


148 


-0 


499 


+0 


364 


-0 


028 


-0 


520 


-0 


182 


+0 


241 


+0 


177 


-0 


348 


-0 


844 


+0 


116 


+0 


397 


-0 


738 


-0 


968 


NGC 


5236: 


XUV 


13 


+0 


280 


+0 


320 


-0 


234 


-0 


554 


+0 


209 


+0 


249 


-0 


202 


-0 


566 


-0 


138 


+0 


201 


-0 


590 


-0 


494 


+0 


241 


+0 


178 


-0 


348 


-0 


843 


+0 


116 


+0 


397 


-0 


738 


-0 


968 


NGC 


5236: 


XUV 


14 


+0 


160 


-0 


130 


-0 


004 


-0 


284 


+0 


261 


-0 


084 


+0 


060 


-0 


237 


+0 


018 


-0 


390 


-0 


347 


-0 


387 


+0 


223 


-0 


277 


-0 


186 


-0 


699 


+0 


192 


-0 


771 


-0 


058 


-1 


202 


NGC 


5236: 


XUV 


17 


+0 


630 


+0 


050 


-0 


394 


-0 


194 


+0 


307 


+0 


031 


-0 


230 


-0 


156 


+0 


364 


-0 


028 


-0 


520 


-0 


182 


+0 


241 


+0 


177 


-0 


348 


-0 


844 


+0 


147 


+0 


077 


-0 


865 


-0 


918 


NGC 


5236: 


XUV 


20 


-0 


050 


-0 


540 


+0 


076 






+0 


104 


-0 


535 


+0 


032 






-0 


201 


-0 


478 


-0 


309 






+0 


154 


-0 


627 


-0 


135 






+0 


192 


-0 


771 


-0 


058 






NGC 


5236: 


XUV 


21 


-0 


300 


-0 


590 


+0 


076 


-0 


474 


-0 


092 


-0 


616 


-0 


018 


-0 


418 


-0 


204 


-0 


679 


-0 


300 


-0 


457 


+0 


154 


-0 


627 


-0 


135 


-0 


801 


+0 


192 


-0 


771 


-0 


058 


-1 


202 


NGC 




XUV 




-0 


290 




270 


+0 


076 


-0 






074 


-0 


426 








440 


— 


194 




456 




310 




468 


+0 


102 


-0 




-0 










192 




771 






-1 


202 



NGC 4625: 
NGC 4625: 
NGC 4625: 
NGC 4625: 
NGC 4625: 
NGC 4625: 
NGC 4625: 



XUV 01 
XUV 02 
XUV 05 
XUV 06 
XUV 07 
XUV 09 
XUV 12 



+0.080 
+0.730 
+0.470 
+0.580 
+0.530 
+0.710 
+0.620 



+0.590 
+0.300 
+0.220 
-0.070 
+0.440 
+0.960 
+0.070 



-0.804 
-0.594 
-0.354 
-0.354 
-0.634 
-0.394 
-0.644 



-0.484 
-0.333 
-0.203 
-0.583 

-0.234 



-0.096 
+0.128 
+0.240 
+0.307 
+0.119 
+0.162 
+0.157 



+0.506 
+0.271 
+0.239 
+0.031 
+0.406 
+0.814 
+0.116 



-0.787 
-0.522 
-0.336 
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-0.483 
-0.385 
-0.468 



-0.372 
-0.312 
-0.156 
-0.510 

-0.410 



+0.141 
+0.386 
+0.364 
+0.364 
+0.361 
+0.722 
+0.364 



+0.547 
+0.333 
-0.028 
-0.028 
+0.528 
+0.888 
-0.028 



-0.802 
-0.709 
-0.520 
-0.520 
-0.754 
-0.337 
-0.520 



H0.141 
-0.182 
-0.182 
-0.100 

-0.182 



+0.240 
+0.204 
+0.241 
+0.277 
+0.204 
+0.223 
+0.105 



+0.654 
+0.387 
+0.178 
-0.146 
+0.387 
+0.918 
+0.038 



-0.806 
-0.751 
-0.348 
-0.217 
-0.751 
-0.783 
-0.690 



-0.954 
-0.843 
-0.838 
-0.954 

-0.947 



+0.197 
+0.116 
+0.116 
-0.041 
+0.116 
+0.221 
+0.147 



+0.690 
+0.397 
+0.397 
-0.027 
+0.397 
+0.917 
+0.077 



-0.748 
-0.738 
-0.738 
-0.739 
-0.738 
-0.735 
-0.865 



-0.968 ' 
-0.968^- 
-1.013*^ 
-0.968 



Note. — (1) Region identification, (2) logarithm of the mcasur ed [OII]AA3726,3729AA/H/3 1 ine ratio (corrected for extinction), (3) logarithm of the measured [OIII] A5007A/H/3 line ratio, (4) 
logarithm of the measured [Nil] A6584A/H/3 line ratio, (5) logarithm of the measured [SII]AA6717,673lAA/H/3 1 ine ratio, (6-9) the same as (2-5) for the line ratios predicted by the best-fitting 
photoionization model assuming a solar N/O ratio and a single star as ionization source, (10-13) the same as (6-9) for models with O/H-scalcd N/O ratio, (14-17) the same as (6-9) for models with 
solar N/O ratio but a starburst as ionization source, (18-21) the same as (6-9) for models with O/H-scalc N/O ratio and a starburst as ionization source. 
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Fig. 1. — False-color RGB image of the southern region of the XUV disk of M83. This 
panel is a composition of the GALEX FUV (blue), ground-based continuum-subtracted Ha 
(green), and ?7-band (red) images. The Ha-selected regions whose spectroscopy we present 
in this paper appear as green compact sources in this figure and are located right below the 
corresponding region identification number (see Table 2). 
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Fig. 2. — False-color RGB image of the XUV disk of NGC 4625. In this case we show a 
combination of the smoothed asinh-scaled FUV image (see Gil de Paz et al. 2005a) (blue), the 
ground-based continuum-subtracted Ha image (green), and the S-band image (red). The 
Ha-selected regions whose spectroscopy we present in this paper appear as green compact 
sources in this figure and are located right below the corresponding region identification 
number (see Table 2). 
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Fig. 3. — Individual spectra of Ha-selected regions in the XUV disk of M83, plus three 
regions located in its optical disk (NGC 5236: XUV 20, 21, and 22). Only regions with 
detected Ha emission at the approximate redshift of M 83 are shown. Identification number 
(according to Table 2) and scaling factors are shown. 
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Fig. 4. — Individual spectra of Ha-selected regions in the XUV disk of NGC 4625. Only 
regions with detected Ha emission at the approximate redshift of NGC 4625 are shown. 
Identification number (according to Table 2) and scaling factors are shown. Spectral regions 
with high sky-subtraction residuals in the spectra of NGC 4625: XUV 06 and NGC 4625: 
XUV 11 have been blanked out for the sake of clarity. 
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Fig. 5. — a) Flux-calibrated sky (+object) spectra within each of the slitlets obtained as 
part of the spectroscopic observations of the XUV disk of M 83 (different colors are used for 
clarity). The absolute differences between the individual spectra are of the order of 0.3 mag. 
b) Differences after small global offsets (the same at all wavelengths but different for each 
spectrum) are applied. The r.m.s. of the individual spectra are now smaller than 0.2 mag at 
almost all wavelengths. c,d) The same as a,b for the observations of NGC 4625. 
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Fig. 6. — Diagnostic diagrams for the emission-line regions in the XUV disks of M83 and 
NGC 4625. a) [OIII]A5007A/H/? versus [NII]A6584A/Ha extinction-corrected line ratios. 
The data for M83 are coded based on their galactocentric distance: regions in inner (optical) 
disk (r<5.5'; filled circles), regions in the inner XUV disk (5.5'<r<10'; filled triangles), and 
regions in the outer XUV disk (r>10'; filled stars). For comparison we show the location 
of the UCM-Survey galaxies, a complete sample of local star-forming galaxies (SFG) and 
AGN selected by their emission in Ha (Gallego et al. 1996), coded by spectroscopic type. 
The lines are the predictions of CLOUDY single-star photoionization models for an age of 
the star of 1 Myr and a solar N/O abundance ratio (see Section [3TB1 . Solid lines correspond 
to models with fixed mass for the ionizing star and metallicities ranging between Z /1O and 
Z while dot-dashed lines correspond to models with fixed metallicity and different masses 
between 20 and 85 M Q . Typical errors on the line ratios are shown at the top right corner 
of the diagram, b) [OIII]A5007A/H/3 versus [OII]AA3726,3729AA/[OIII]A5007A extinction- 
corrected line ratios. Here we also show the predictions of photoionization models for evolving 
starbursts (Stasinska & Leitherer 1996), with different lines showing the predictions for 
different metallicities (see text for more details on the models), c) The same as a for 
NGC 4625. d) The same as b for NGC 4625. The starburst photoionization models shown 
here correspond to one-tenth-solar abundance models with two different densities, 10 and 
10 4 cm~ 3 . 
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Fig. 7. — a) Radial distribution of the nebular oxygen abundance in M83. Filled dots 
correspond to metallicities derived from the line ratios measured in this work. Open circles 
are the metallicities derived using the line ratios measured by published by Webster & 
Smith (1983) in six HII regions in the inner optical disk of M83. The black (grey) dots 
are the metallicities obtained assuming the low (high) metallicity branch for the calibration 
of the R23 parameter for regions NGC 5236: XUV 09 to 14. The black (grey) dashed- 
line is the corresponding best linear fit. The solid-line is the best fit obtained after using 
the metal abundance information provided by the best-fitting single-star photoionization 
models (open squares; see text for details). These oxygen abundances were derived using 
the McGaugh (1991) calibration of the R23 parameter. Coefficients for all these fits and 
those obtained using the calibration of Pilyugin & Thuan (2005) are given in the text, b) 
Radial distribution of the internal extinction in H/5 [C(H/3)] for M83. We have adopted the 
Galactic color excesses in the direction of M83 and NGC 4625 to be 0.066 and 0.018 mag, 
respectively (Schlegel, Finkbeiner, & Davis 1998). Symbols have the same meaning as in 
panel a. c,d) The same as a,b for NGC 4625. 



